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We develop a new ab initio many-body approach capable of describing simultaneously both bound 
and scattering states in light nuclei, by combining the resonating-group method with the use of 
realistic interactions, and a microscopic and consistent description of the nucleon clusters. This 
approach preserves translational symmetry and Pauli principle. We outline technical details and 
present phase shift results for neutron scattering on ^H, ''He and ^"Be and proton scattering on 
'''*He, using realistic nucleon-nucleon (NN) potentials. Our A — i scattering results are compared 
to earlier ab initio calculations. We find that the CD-Bonn A'^A^ potential in particular provides an 
excellent description of nucleon-*H[e S-wave phase shifts. On the contrary, the experimental nucleon- 
■*He P-wave phase shifts are not well reproduced by any A'^A'^ potential we use. We demonstrate 
that a proper treatment of the coupling to the n-^°Be continuum is successful in explaining the 
parity-inverted ground state in ^^Be. 



I. INTRODUCTION 

Nuclei are open quantum systems with bound states, 
unbound resonances, and scattering states. A realistic 
ab initio description of light nuclei with predictive power 
must have the capability to describe all the above classes 
of states within a unified framework. Over the past 
decade, significant progress has been made in our un- 
derstanding of the properties of the bound states of light 
nuclei starting from realistic nucleon-nucleon (NN) inter- 
actions, see e.g. Ref. [Ij and references therein, and more 
recently also from A'' A'' plus three-nucleon (NNN) inter- 
actions [H [21 H] . The solution of the nuclear many-body 
problem is even more complex when scattering or nu- 
clear reactions are considered. For A = 3 and 4 nucleon 
systems, the Faddeev [5] and Faddeev-Yakubovsky [6] 
as well as the hyperspherical harmonics (HH) [7\ or the 
Alt, Grassberger and Sandhas (AGS) [8^ methods are ap- 
plicable and successful. However, ab initio calculations 
for scattering processes involving more than four nucle- 
ons overall are challenging and still a rare exception [UJ. 
The development of an ab initio theory of low-energy nu- 
clear reactions on light nuclei is key to further refining 
our understanding of the fundamental nuclear interac- 
tions among the constituent nucleons and providing, at 
the same time, accurate predictions of crucial reaction 
rates for nuclear astrophysics. 

Recently we combined the resonating-group method 
(RGM) [TOl El m [SI m IS and the ab initio no- 
core shell model (NCSM) [TB], into a new many-body 
approach [T7j {ab initio NCSM/RGM) capable of treat- 
ing bound and scattering states of light nuclei in a unified 
formalism, starting from the fundamental inter-nucleon 
interactions. The RGM is a microscopic cluster tech- 
nique based on the use of A-nucleon Hamiltonians, with 
fully anti-symmetric many-body wave functions built as- 
suming that the nucleons are grouped into clusters. The 
NCSM is an ab initio approach to the microscopic cal- 
culation of ground and low-lying excited states of light 
nuclei with realistic two- and, in general, three-nucleon 
forces. The use of the harmonic oscillator (HO) ba- 



sis in the NCSM results in an incorrect description of 
long-range correlations and a lack of coupling to contin- 
uum. The first NCSM applications to nuclear reactions 
required a phenomenological correction of the asymptotic 
behavior of overlap functions [18^ . The present approach 
is fully ab initio. We complement the ability of the RGM 
to deal with scattering and reactions with the use of re- 
alistic interactions, and a consistent ab initio description 
of the nucleon clusters, achieved via the NCSM. Within 
this new approach we studied the n-'^H, n-'*He, n-^°Be, 
and p-'^'^He scattering processes, and addressed the par- 
ity inversion of the ^^Be ground state (g.s.), using real- 
istic A^A^ potentials. In this paper, we give the technical 
details of these calculations, discuss results published in 
Ref. [IT more extensively and present additional results. 

In Scct.|ll] we present technical details of our approach. 
We give two independent derivations of the NCSM/RGM 
kernels, we discuss orthogonalization of the RGM equa- 
tions and give illustrative examples of the kernels. Re- 
sults of ab initio NCSM/RGM applications to A = 4, 
A = 5 and A= 11 systems are given in Sect. |III[ Conclu- 
sions are drawn in Sect. |IV| and some of the most complex 
derivations are summarized in Appendix 



II. FORMALISM 

The wave function for a scattering process involving 
pairs of nuclei can be cast in the form 



I* 



E 



(1) 



through an expansion over binary-cluster channel-states 
of total angular momentum J, parity tt, and isospin T, 
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The internal wave functions of the colliding nuclei (which 
we will often refer to as clusters), contain A— a and a 
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nucleons {a<A), respectively, are antisymmetric under 
exchange of internal nucleons, and depend on transla- 
tionally invariant internal coordinates. They are eigen- 
states of and if(a)i the (A~a)- and a-nucleon 

intrinsic Hamiltonians, respectively, with angular mo- 
mentum quantum numbers Ii and I2 coupled together 
to form channel spin s. For their parity, isospin and 
additional quantum numbers we use, respectively, the 
notations T^i,Ti, and a^, with i = 1,2. The channel 
states ([2| have relative angular momentum t. Denoting 
with {ri,i — 1,2,-- - ,A} the A single-particle coordi- 
nates, the clusters centers of mass are separated by the 
relative vector 



'^A—a,a — ^A—a,a^A — a.a 
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The symbols Yg and 6 denote a spherical har- 
monic and a Dirac delta, respectively. The inter- 
cluster anti-symmetrizer for the {A— a, a) partition 
m Eq. M can be schematically written as Ai, = 
aTjAl]^/^' 



[{A-a)\ 



'^p(— )^'P, where P are permutations 



among nucleons pertaining to different clusters, and p 
the number of interchanges characterizing them. 

The coefficients of the expansion with respect to the 
channel index i/ = {j4—oai/]"^Ti; aa2l2^T2\ s£} are the 
relative-motion wave functions g'l '^{r), which represent 
the unknowns of the problem. They can be determined 
by solving the many-body Schrodinger equation in the 
Hilbert space spanned by the basis states A,^ \^'^ 



,.rT\ 



^ j drr^ 
where 
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Au'HAu 



= 0, 
(4) 

(5) 
(6) 



are called the Hamiltonian and norm kernels, respec- 
tively. Here E is the total energy in the center-of-mass 
(cm.) frame, and H is the intrinsic A-nucleon micro- 
scopic Hamiltonian, for which it is useful to use the de- 
composition, e.g.: 



H ^T, 



\el{r) + Vrel + Vc{r) + Hf^^-a) + H^a) ■ (7) 



Further, Trei(r) is the relative kinetic energy and Vrci 
is the sum of all interactions between nucleons belong- 
ing to different clusters after subtraction of the average 
Coulomb interaction between them, explicitly singled out 
in the term Vc(^) = Zii,Z2i^e'^ /r, where Zi^ and Z2u are 



the charge numbers of the clusters in channel v. 
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In the above expression we explicitly distinguished 
between nucleon-nucleon, nuclear (Vat) plus Coulomb 
(point and average), and three-nucleon (V^^^^^^) com- 
ponents of the inter-cluster interaction. The contribution 
due to the nuclear interaction vanishes exponentially for 
increasing distances between particles. Thanks to the 
subtraction of Vc{r), the overall Coulomb contribution 
presents a r^^ behavior, as the distance r between the 
two clusters increases. Therefore, Vrei is localized also 
in presence of the Coulomb force. In the present paper 
we will consider only the NN part of the inter-cluster 
interaction, and disregard, for the time being, the term 
V?{^ X. The inclusion of the three-nucleon force into 

(A— a, a) 

the formalism, although more involved, is straightfor- 
ward and will be the matter of future investigations. Fi- 
nally, although in Eq. (|8| the strong part of the NN 
force (Vat) is represented as a local potential, the above 
separation of the Hamiltonian as well as the rest of the 
formalism presented throughout this paper are valid also 
in the presence of a non-local potential. 



A. Cluster eigenstate calculation 

We obtain the cluster eigenstates entering Eq. ^ 
by diagonalizing H(^A-a) and in the model space 

spanned by the NCSM basis. This is a complete HO 
basis, the size of which is defined by the maximum num- 
ber, A^max, of HO quanta above the lowest configuration 
shared by the nucleons (the definition of the model-space 
size coincides for eigenstates of the same parity, differs 
by one unity for eigenstates of opposite parity; the same 
HO frequency is used for both clusters). If the NN (or 
NNN) potential used in the calculation generates strong 
short-range correlations, which is typical for standard ac- 
curate NN potentials, the and H{a) Hamiltoni- 
ans are treated as NCSM effective Hamiltonians, tailored 
to the A'max truncation, obtained employing the usual 
NCSM effective interaction techniques [4, 16J. The ef- 
fective interactions are derived from the underlying NN 
and, in general, three-nucleon potential models (not in- 
cluded in the present investigations) through a unitary 
transformation in a way that guarantees convergence to 
the exact solution as the model-space size increases. On 
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the other hand, if low-momentum NN potentials, which 
have high-momentum components already transformed 
away by unitary transformations, are employed in the cal- 
culations, the and ^^(a) Hamiltonians are taken 
unrenormalized or "bare." 

Thanks to the unique properties of the HO basis, we 
can make use of Jacobi-coordinate wave functions [151 HO] 
for both nuclei or only for the lightest of the pair (typi- 
cally a < 4) referenced further on as projectile, and still 
preserve the translational invariance of the problem. In 
the second case we expand the eigenstates of the heav- 
ier cluster (target) on a Slater-determinant (SD) basis, 
and remove completely the spurious cm. components in 
a similar fashion as in Refs. [HI |2TJ |22]. We exploited 
this dual approach to verify our results. The use of the 
SD basis is computationally advantageous and allows us 
to explore reactions involving p-shell nuclei. 



where x — \/ ^{ri ~ and p = i/ ^(pi — P2), we intro- 



2V' i ■ ^/ 1^"^ ^ ~ Y 2' 
duce here a second, modified two-nucleon Hamiltonian, 
deprived of the nuclear interaction: 



— tlQ2 



V{2 = 



(10) 



1 

2m ' 2' A 
The modified two-body effective interaction is then deter- 
mined from the two-nucleon Hermitian effective Hamil- 
tonians i?2off and -ff2cff ' obtained via the Lee-Suzuki sim- 
ilarity transformation method [23] starting from Eqs. ^ 
and (10), respectively: 



^2off — H2cS - H2 



U- (11) 

We note that i) Vj'^fj Vn in the limit N^mx ^ 00, and 
ii) for each model space, the renormalizations related to 
the kinetic energy and the HO potential introduced in 
H2cS are compensated by the subtraction of -ff2eff- 



B. 



Interaction between nucleons belonging to 
different clusters 



In calculating (5j6 1, all "direct" terms arising from the 
identical permutations in both Ai, and A^> are treated 
exactly (with respect to the separation r) with the excep- 
tion of {^'l^^l Vrei \ The latter and all remaining 
terms are localized and can be obtained by expanding 
the Dirac 6 of Eq. ([2]) on a set of HO radial wave func- 
tions with identical frequency fl, and model-space size 
-^max consistent with those used for the two clusters. 
The rate of convergence of these terms is closely related 
to the nuclear force model adopted in the Hamiltonian 
([7]). For most nuclear interaction models that generate 
strong short-range nucleon-nucleon correlations the large 
but finite model spaces computationally achievable are 
not sufficient to reach the full convergence through a 
"bare" calculation. In these cases it is crucial to uti- 
lize effective interactions tailored to the truncated model 
spaces. In our approach the effective interactions are de- 
rived from the underlying NN potential through a uni- 
tary transformation as already pointed out in the pre- 
vious Subsection. While the cluster eigenstates are ob- 
tained employing the usual NCSM effective interaction 
[TS] , in place of the bare NN nuclear potential Vn en- 
tering Vroi ([8| we adopt a modified two-body effective 
interaction, V2^g, which avoids renormalizations related 
to the kinetic energy. While the kinetic-energy renormal- 
izations are appropriate within the standard NCSM, they 
would compromise scattering results obtained within the 
NCSM/RGM approach, in which the relative kinetic en- 
ergy and the average Coulomb interaction between the 
clusters are treated exactly. More specifically, in addi- 
tion to the relevant two-nucleon Hamiltonian (see also 
Refs. [Ulin]) 
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C. Coordinates and basis states 

We neglect the difference between proton and neu- 
tron masses, and denote the average nucleon mass with 
m. The formalism presented in this paper is based 
both on the single-particle Cartesian coordinates, {r^, i — 
1, 2, • • • , A}, and on the following set of Jacobi coordi- 
nates: 

^o-JIj:^^ (12) 

the vector proportional to the center of mass (cm.) co- 
ordinate of the A-nucleon system (i?c.m. — :^$o); 



= \/ 2(^1 -^2) 



k + 1 



1 



i=l 



,2 < fc < A-a-1 ; 



(13) 

the translationally-invariant internal coordinates for the 
first A^a nucleons; 



VA-a = 



(A - a)a 
A 



A— a 



A- a 



i=A-a+l 



(14) 

the vector proportional to the relative position between 



the cm. of the two clusters (f^-a,a = y {A-a)a VA-a); 

and, finally, 

"1 

V E ^-4-^+1 - ^A-k 



'A-k 



^A- 



k+1 
1 



i=l 



,a-l > fc > 2, 



(15) 



the translationally-invariant internal coordinates for the 
last a nucleons. 
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1. Jacobi basis 

Nuclei are translationally invariant systems. There- 
fore, the use of Jacobi coordinates and translationally- 
invariant basis states represents a "natural" choice for 
the solution of the many-nucleon problem. 

Working with the Jacobi relative coordinates of 
Eqs. (13 1, (14 1, and (15 1, it is convenient to introduce 



the (translationally-invariant) Jacobi channel states 



\A-aaiI,^'T,)\aa2lpT2)) 

(.rX) S{j^ - rjA-a) 



(sT) 



VVA-a 



(16) 



which are clearly proportional to the binary-cluster basis 
presented in Eq. ([2]): 



I1> 



{A-a)a 
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l<i> 
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The clusters intrinsic wave functions depend on their re- 
spective set of Jacobi, spin (tr) and isospin (r) coordi- 
nates: 

(6- • • iA-a-l(JV ■ ■ CTA-aTl- ■ ■ TA-a\A-a ail'^^Ti) , (18) 



■TA\aa2^'T2) 



(19) 

and are obtained by diagonalizing the H(^A-a) find i?(a) 
intrinsic Hamiltonians in the model spaces spanned by 
the NCSM Jacobi-coordinate basis |20j. The same HO 
frequency is used for both clusters. The model-space 
size coincides for eigenstates of the same parity and dif- 
fers by one unit for eigenstates of opposite parity. 

In calculating the integral kernels of Eqs. ([s]) and ([6|, 
{^■l,^\Vrei\^ir'^) and all "exchange" terms, arising 
from the permutations in Ai, or A,yi different from the 
identity, are obtained by expanding the Dirac S of Eq. ^ 
on a set of HO radial wave functions with identical fre- 
quency f2, and model-space size iVmax consistent with 
those used for the two clusters: 



l<i> 



{A-a)a 



-1 3/2 



A 



\Rni{r, b) 1$ 



(20) 
(21) 



where the HO Jacobi-channel states are given by 



1$ 



\A-aaiI^'Ti) \aa2IpT2) 

Yl{f]A-a) Rnl{rA-a,a,b) (22) 
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Note that the HO basis states depending on the Ja- 
cobi coordinates introduced in Section III CI are all char- 
acterized by the same oscillator-length parameter 60 = 
y/h/mil. However, the oscillator-length parameter asso- 
ciated with the separation r between the centers of mass 
of target and projectile is defined in terms of the reduced 
mass /i = [(A—a)am\/A of the channel under considera- 
tion: b — y/h/ fifl — y/ A/[(A~a)a\ bo. In the following we 
will drop the explicit reference to the HO length param- 
eter in the arguments of the HO radial wave functions, 
and in the HO Jacobi channel states |<i>;!„-^). 



2. Single-particle Slater- determinant basis 

Thanks to the unique properties of the HO basis, we 
can make use of Jacobi-coordinate wave functions ^f^EU] 
for both nuclei or only for the lighter of the pair (typically 

^ 4), and still preserve the translational invariance of 
the problem (see also discussions in Refs. 21, 22J. In 
the second case we introduce the SD channel states 
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in which the eigenstates of the (A— a)-nucleon fragment 
are obtained in the SD basis. 



(fi- • -rA-aCTr --(JA-aTv ■ -TA-alA-aaill^Ti] 



SD, 



(25) 



i.e., by using a shell-model code (such as e.g. Antoine jH] 
or MFD [25]), and contain therefore the spurious motion 
of the (A— a)-nucleon cluster cm. The SD and Jacobi- 
coordinate eigenstates are related by the expression: 

|A-aai/iri)sD = \A-aaJ^T^) ^oo(i?ii:"') • (26) 
The cm. coordinates introduced in Eqs. ([24|) and (p6| 



1 



- A — a 
i=l 
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(27) 

are an orthogonal transformation of the c m. a nd relative 
coordinates of the A-nucleon system, ^0 (12 1 and ffA-a 
( 14 1, respectively: 



VA-a = 



Co 




A-c 
~~A 



]^{A-a) 



(28) 



(29) 



Therefore, in the SD basis of Eq. (24 1, the HO wave func 



tions depending on these coordinates transform accord- 
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ing to 



1 2 ■ • • A-2 A 



A-1 



J2 {QQnU\n4rNU)^^{^nj,XVA-a) VNL{ia)y ^ , 

(30) 

where the coefficients of the expansion are generalized 
HO brackets for two particles with mass ratio d = 
that can be calculated as described e.g. in Ref. [35]. As 
a result the SD and Jacobi channel states are related by: 
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FIG. 1: (Color online.) Diagrammatic representation of the 
"direct" (a) and "exchange" (fo) components of the norm ker- 
nel. The first group of circled black lines represents the first 
cluster, the bound state of A—1 nucleons. The separate red 
line represents the second cluster, in the specific case a single 
nucleon. Bottom and upper part of the diagram represent 
initial and final states, respectively. 



where = {A—aaiIiTi\ a 0:2/2 ?2; s£r} ■ It is therefore 
possible to extract the translationally-invariant matrix 
elements from those calculated in the SD basis, which 
contain the spurious cm. motion, by inverting the fol- 
lowing expression: 



In calculating (|5| and ([6|, it is convenient to isolate 
the "direct" terms arising from the identical permuta- 
tion in A. Considering that the full A-nucleon Hamil- 
tonian commutes with the inter-cluster anti-symmetrizer 
{[A,H] = 0), and that 



sD(*;;rj Ot.i 



SD 



A-1 



A^ 



(34) 
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X {nr£rNLe\OOnU)^ {n'XNL£\00n'£'£'' 



(32) 



where Ot.i. is any scalar and parity-conserving 
translational-invariant operator (Ot.i. = A, AHA, etc.). 

We exploited this dual approach to verify our results. 
The use of the SD basis is computationally advantageous 
and allows us to explore reactions involving p-shell nuclei. 



we can write the following expression for the norm kernel 
in the SNP basis: 

K'Jir', r) = 5,, . ^^^^ + K%{t\ r) . (35) 

Here, we have singled out the non-local exchange part of 
the matrix elements in the term (we drop for simplicity 
the J'^T superscript) 



A-1 



(36) 



= -{A - l)Y,Rn'f{r')Rni{r) 



D. Translational invariant kernels in the 
single-nucleon-projectile basis 

All calculations in the present paper were carried out 
in the single-nucleon projectile (SNP) basis, i.e., using 
binary-cluster channels ([2| with a—1. In this case, the 
■d coordinates are not defined, the channel index reduces 
to ^ = {A—1 ail^^Ti] s£}, and the inter-cluster 

anti-symmetrizer is simply given by 



V3 



A-1 



(33) 



X (^iZ\pA-l^A 



) . (37) 



In deriving Eq. (37) we used the expansion (21), and 



took advantage of the internal symmetry properties of the 
(A— l)-cluster wave function. A similar decomposition 
can be performed also for the Hamiltonian kernel. 



A-1 



^ i/'r' 



H 



1=1 



(38) 



(39) 
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FIG. 2: (Color online.) 
Diagrammatic repre- 
sentation of "direct" (c 
and d) and "exchange" 
(e) components of the 
potential kernel (see 
also Caption of Fig. [T]). 



where we divided {^■lJJ\VrciA'^ \^ir'^) into "direct" 
and "exchange" potential kernels according to: 

V°,(/,r) = {A-l)Y,Rn'i'{r')Rni{r) 



x{<^i:;^\VA-i,Aa-PA-K 



(40) 



V.%{r\r) = -{A-l){A-2)J2Rn'i'{r')Rni{r) 



X ( ^-TT 



Pa-i,a Va-2,A-1 



(41) 



As pointed out in Sec. |TT] the channel states ([2| are not 
anti-symmetric with respect to the exchange of nucle- 
ons pertaining to different clusters (fully anti-symmetric 
states are recovered through the action of the opera- 
tor A^). As a consequence, the Hamiltonian kernel as 
defined in Eq. (38 1 is explicitly non Hermitian. Using 



AHA ~ \{A^H + HA?), we introduce the Hermitized 
Hamiltonian kernel HLJ in the form 



A-l 



ni:j{r',r)=Ui:^ h-\y. {HhA+hAH) 



(42) 

Finally, we note that, according to Eqs. (7[8) and 
Eqs. ( 35|37 l, the contribution of the average Coulomb po- 
tential to the Hermitian Hamiltonian kernel ( 42 ) amounts 
overall to: 



2'5.v[Vb(r') + VcM] 



5{r' - r) 



(43) 



1. J acohi- coordinate derivation 



The main technical as well as computational challenge 
of the NCSM/RGM approach lies in the evaluation of 
norm and Hamiltonian kernels. The analytical expres- 
sions for the integral kernels of Eqs. (37), (40), and (41) 



assume a particularly involved aspect in the model space 
spanned by the HO Jacobi channel states of Eq. ( [22| . 
Here we discuss the exchange-part of the norm kernel for 
the ^ = 3 system (a = 1), which is representative of the 



Jacobi-coordinate formalism without requiring overly te- 
dious manipulations. Interested readers can find a com- 
pilation of all Jacobi-coordinate formulae, along with an 
outline of their derivation, in Appendix \K\ 

The HO Jacobi channel state of Eq. (|22]) for the (2, 1) 
partition can be written as 



1$ 



1 Is 



X I [(m^iSi/iTi; --)sT; n£] J^T) , (44) 

where we have expanded the two-nucleon target wave 
function onto HO basis states depending on the Jacobi 
coordinate defined in Eq. ( 13 1 



(^1 (Tl 0-2 Tl r2 1 71. 1 €l S 1 /l Tl ) 



(45) 



and (ni^iSi/iTi|2ai/fiTi) are the coefhcients of the ex- 
pansion. Here ni,ii 'Ave the HO quantum numbers cor- 
responding to the harmonic oscillator associated with ^i, 
while si,/i, and Ti are the spin, total angular momen- 
tum, and isospin of the two-nucleon channel formed by 



nucleons 1 and 2, respectively. Note that the basis (45) is 



anti-symmetric with respect to the exchange of the two 
nucleons, (-)'^i+*i+Ti = -l. 



According to Eq. (37), in order to obtain the exchange 



part of the norm kernel we need to evaluate matrix ele- 
ments of the permutation corresponding to the exchange 
of the last two particles, in this case P23- This task can 
be accomplished by, e.g., switching to a more convenient 
coupling of the three-nucleon quantum numbers 

|[(nifiSi/iri;^^)sT;n^]j"r 
= ^Z/i(-)^^+^^ + 5+J^^i I 

Z L 2 J 

x5:A^(-)-+^+^+4f 5 1 

A 

X I [(m^i , n£)K; {s, ^)Z] | {T, ^)t) , (46) 

and observing that, as a result of the action of 
P23, the HO state {^iff2\{ni£i,n£)A) is changed into 
{£^lf]2\{ni£i,n£)A). The new set of Jacobi coordinates 
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£^1 and (obtained from ^1 and 772, respectively, by ex- ones. Consequently, the HO states depending on them 
changing the single-nucleon indexes 2 and 3) can be ex- are related by the orthogonal transformation 
pressed as an orthogonal transformation of the unprimed 



(^i'%'IK4,«^)A)= J2 (-)^+^^"^(A^^,AAi£i,A|ni£i,n^,A)3(eir/2|(AAiA,iVL)A), (47) 

where the elements of the transformation are the general HO brackets for two particles with mass ratio d — 3. 

After taking care of the action of P23 also on the spin and isospin states, one can complete the derivation and write 
the following expression for the A = 3 exchange part of the norm kernel in the SNP basis: 



e[ z s' ^ ( ii z s ) ( £1 z s 



n'e',n[£[,A\niii,ni,A)3. (48) 



Here we remind that the index v stands for the collection 
of quantum numbers {A—1 ail^'^Ti; l^; s£}, while v' 
is an analogous index containing the primed quantum 
numbers. 

The derivation of "direct"- and "exchange" -potential 
kernels, although complicated by the need for additional 
orthogonal transformations and the presence of the two- 
body matrix elements of the interaction, proceeds along 
the same lines presented here (see Appendix All. As 
final remark, we note that while the exchange part of 
the norm kernel ( 48 ) and the direct potential kernel ( A4 ) 



are symmetric under exchange of prime and unprimed in- 
dexes, and primed and unprimed coordinates, the same is 



not true of the exchange part of the potential kernel ( A5 1 
Indeed, as anticipated in Sec. HD| the Hamiltonian ker- 
nel defined in Eq. (38 1 is explicitly non Hermitian. 



2. Smgle-particle Slater determinant derivation 



The matrix elements of the operators Pa-i,Aj 

Va-i,a(1 - Pa-i,a), and Pa~i.aVa-2,a-i can be more 
intuitively derived working within the SD basis of 
Eq. ( 24 ) . Using the second-quantization formalism, they 



can be related to linear combinations of matrix elements 
of creation and annihilation operators between (A—l)- 
nucleons SD states. These quantities can be easily cal- 
culated by shell model codes. Here we outline the main 
stages of the derivation. 



The SD basis (24) simplifies in the case of a single- 



nucleon projectile to 



l*;^7)sD 



[( 



\A-laihTi] 



SD 



ll\)r^Y,irA)f''R..irA) 



I J J 



1^-1 "1-^1^1)30 ^neairACTATA 



1 (J^T) 



(49) 



r 



with u = {^-lai/f^Ti; 1 ii; s£} and the HO 
single-particle wave function y^nejm^^mti^AO'ATA) — 



Rne{rA){yeifA)xi{'^A))^;l\^mA'rA)- To obtain the 



exchange part of the norm kernel ( 37 1 we first 



calculate the permutation operator matrix elements 
within the basis ( 49 1 . By expressing the position 
state of the nucleon {A— I) as I'ta-iO'a-ita-i) = 



{rA-l(JA-lTA- 



we ar- 



8 



rive at 



SD«3|PA,A-I|<i>f/)SD 



X sD(A-la7;r{|||(a|^,^.,a„,,,^.,.)(^-)|||A-la/iri)sD. (50) 

Here, sD(^-la'/ir{ | ||(a^ , . i a„,^,,-, i )(^^' 1 1| A-la/iTi)sD are one-body density matrix elements (OBDME) of the 
target nucleus and a„,^/j/„/ i„j/ = (—1)-' ^"^ +5^™* o„/£,j/_„j/ . Next we extract the corresponding translationally- 



invariant matrix elements, (*i>Jy^„/^'^'''^'' ^1^4.^1-1 1'^lvn^'^''"''^ by inverting Eq. (32l for a— \ and Ot.i. = Pa-\.a- 



The final step follows easily from Eq. ( 37 1 



The same procedure is applied also for calculating "direct"- and "exchange" -potential kernels. In this case the 
transition matrix elements on the SD basis are respectively: 



X 



f J J [f J Jo [ kT 



1 T i 



11^ 



(nij) 



X {{nJJa \ ) {n'i'f i ) JoTo \V\{n£j ^ ) {nthjt ^ ) JqTo) 



and 



sd($;^'^|-Pa,a-iVa-2,a-i|<&;^„^)sd 
1 



E E E E E E Ss'jj'KTkafaKcdTcd 



2(A-l)(A-2) ^ 

jj'Kr rialaja rLbhjb nJcjc ridldjd KaTaK^aTcd 

y^(^^lY'l+]' + 'l+K+]+Ja+ja+Jdl^_l-yri + l^+T + T 

\h\s\ \l[\s'\ illKTA (KaK,dK] (TiTT[] (vTaT^d 



2^ 2j L22 2 



^ Jj J 1 ^' J/J 1/ J J J 1 / J Jaj 

X + ^(naiaja),("'f'i')\/l + ^("c /c J c ) , ("d /rf J d ) ( ("'^'j' ^ ) (^^o^a Ja ^)KcdTcd \V\{nJdjd ^ ) {nclcjJj)KcdTcd) 

^MA-la'I[Tm{al^,/^^^^^^^^^ (52) 



While the "direct" matrix element (51 1 depends on the mer involves only a single nucleon of the target, while the 



OBDME, the "exchange" matrix element (52 1 depends latter involves two nucleons of the target, see also Fig. |2] 
on two-body density matrix elements (TBDME) of the We note that the two-body matrix elements of the inter- 
target nucleus. This is easily understandable as the for- action V are evaluated using just the first two terms of 



9 



Eq. (81, i.e. Vu 



VN{^3) 



as the aver- 



age Coulomb interaction is taken care of with the help 
We also note that as a consistency check, it 



(431 



of Eq 
is possi 

expression (|51| or the expression (|52| by setting the NN 



jle to recover the expression ( 50 1 from either the 



interaction operator V to identity. 



3. Illustrative examples 

The n-a system provides a convenient ground to ex- 
plore the characteristic features of the integral kernels 
obtained applying the NCSM/RGM approach within the 
SNP formalism. Thanks to the tightly-bound structure 
of ^He, an expansion in n-a channel states allows us to 
describe fairly well the low-energy properties of the ^He 
system. The latter (likewise ^Li) is an unbound system, 
its ground state being a narrow P-wave resonance in the 
I i channel. 

Figures [3] to [?) and Table |l] present results of single- 
channel calculations carried out using n-a cluster chan- 
nels with the a particle in its g.s. (note that through- 
out this Section the index v = {4g.s. 0'*'0; 1^^^; can 
and will be simply replaced by the quantum number 
£). The interaction models adopted are the N^LO NN 
potential [27] derived within chiral effective-field theory 
(xEFT) at the next-to-next-to-next-to-leading order, and 
the V\ov/k NN potential derived from AV18 with cut- 
off A — 2.1 fm~^. Although xEFT forces are known to 
present a relatively soft core, the large but finite model 
spaces computationally achievable are still not sufficient 
to reach a full convergence through a "bare" calculation. 
Therefore, for this potential we utilize two-body effec- 
tive interactions tailored to the truncated model spaces 
as outhned in Sec. |IIB[ Results for the Viowfc potential 
are obtained using the "bare" interaction. 

The overall convergence behavior of the integral ker- 
nels is influenced by both the convergence of the eigen- 
states entering the binary-cluster basis, in the specific 
case the '*He g.s., and the convergence of the radial ex- 
pansion of Eq. (21 1. As an example. Fig. [3] presents the 
behavior of the exchange part of the norm kernel with re- 
spect to the increase of the model-space size obtained for 
the J^T = 5^5, and | ^ five-nucleon channels, using 
the N'^LO potential. The corresponding convergence pat- 
tern for the a-particle g.s. energy is shown in the inset. 
In order to allow for the calculation of both positive- and 
negative-parity five-nucleon channels, for a given trun- 
cation iVinax in the /f^Ti = 0+0 model space used to 



1. 



expand the g.s., a complete calculation of Eq. (37) re 
quires an expansion over n-a J^T states up to N^^ 
This is the origin of the odd A^^ax values in the legend of 
Fig. [3] (and following). As we can see from the figure, the 
HO frequency hH. — 19 MeV enables a quite satisfactory 
convergence of both ^He g.s. and n-a radial expansion, 
and hence of the integral kernel. As an example, for the 




FIG. 3: (Color online.) Dependence on A'max of the "ex- 
change" part of the diagonal norm kernel for the n-''lle(g.s.) 
1^1 (^Si/2), and I I CP-i/2) channels as a function of the 
relative coordinate r at r' = 1 fm, using the N'^LO A'' A'' po- 
tential .27il at hQ — 19 MeV. In the inset, convergence pattern 
of the energy of the ^He g.s., used to build the binary-cluster 
basis. The green dashed line indicates the previous NCSM 
evaluation of Eg.s. = -25.39(1) MeV [29]. 



3% or less off the converged (Af„iax = 21) curve in the 
whole r-range up to 4.5 fm. An analogous analysis of the 
kernels yields a somewhat larger relative difference 
(less than 10%) between N^^-^ = 17 and 21 in the range 
between 1 and 4 fm, while the discrepancy increases to- 
wards the origin. In this regard, we note that the | ^ 



one. 



kernel overall is an order of magnitude smaller than the 

1+1 
2 2 

The convergence rate for V\o 
Fig. 4| is clearly much faster. Here the results 



51/2 channel the N^ 



17 result is already within 



(see upper panel of 

°3/2 

for tlie two largest model spaces (A^max =15 and 17) are 
within 0.5% or less in the whole region up to 5 fm. 

Despite the mild differences in magnitude and strength 
distribution for small r, r' values, the and ^^3/2 re- 

sults of Figs. |3]and |4]present essentially the same shape, 
and same range of about 5 fm. This can be observed also 
in Fig.js] which shows once again the "^81/2 partial wave, 
in terms of contour plots (note that the 5*1/2 curves of 
Figs. |3] and |4] correspond to slices of the current plot 
along the r' = 1 fm line). In particular it is clear that 
the ^51/2 kernels for the two different A^A^ potentials as- 
sume almost-identical values starting from r, r' = 2 fm, 
the N'^LO results being much shallower near origin and 
overall less symmetric than those obtained with Viowfe- 
The latter features reveal differences in the structure of 



10 



0.05 

0.04 

0.03 

"^0.02 

^0.01 

i 

-0.1 
-0.2 
-0.3 
-0.4 
-0.5 
-0.6 



i 1 i 1 i 1 i 


1 ' 1 


1 ' 1 ' 

n + cx{g.s.) ■ 


. / \ f 1 


N 

^ ^ max 
.... 9 

11 

— ■ 13 

15 

17 




.1.1.1. 


-v,^ no. = 18 MeV . 




/' 

/ 




^lowk 




/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ N —17 

y ' max — ' 


1^3 (^-^i/a) 

— fU'Py^) 

— 2 2 y ^3/2) 

2 2 ( -^3/2) 




I.I.I. 


1 1 


1 . 1 









4 

r [fm] 



FIG. 4: (Color online.) "Exchange" part of the diagonal norm 
kernel for the n-'*He(g.s.) and channels as a func- 

tion of the relative coordinate r at r' = 1 fm, using the Viowfe 
NN potential gS] at ftfi = 18 MeV. The upper panel shows 
the model-space dependence of the ^P3/2 component. 



tions of higher relative angular momenta are of the same 
order or smaller than the ^I?3/2 partial wave. It is appar- 
ent that the ^5*1/2 channel dominates all over the others 
and is negative. This is an effect of the Pauli exclu- 
sion principle, which forbids more than four nucleons in 
the s-shell of a nuclear system. The four nucleons form- 
ing the ^He g.s. sit mostly in the Ohil shell. Accord- 
ingly, in the ^6*1/2 channel the "exchange" -part of the 
norm kernel suppresses the (dominant) Ohfl contribution 

to the S function of Eq. (l35l (and, consequently, to the S- 

— 1+1 
wave relative-motion wave function 5|^g^ ) coming from 
the fifth nucleon in s-shell configuration. More precisely, 

the diagonalization of the "exchange" part of the norm 

1 + 1 

kernel reveals the presence of an eigenvector gg p ^ with 
eigenvalue 7r — —1, i.e. a Pauli-forbidden state: 



drAA-(r',r) 



1 + 1 
-2 2 



M-7r0oV'('^') 



(53) 



Table |T]presents the three largest-negative eigenvalues for 
the adopted A'^A^ potentials along with their dependence 
upon the model space size. For both interactions the 
first eigenvalue clearly corresponds to a Pauli-forbidden 
state. Once again, the rate of convergence for Viowk is 
visibly faster than for N^^LO, and, despite the differences 
noted in the integral kernels, the overall results for the 



the a particle obtained within the N^^LO and Viowfc A^A^ 
interactions. We note that g.s. energy and point-proton 
root-mean-square radius of the a particle are —25.39(1) 
MeV, 1.515(2) fm and -27.77(1) MeV, 1.4239(2) fm with 
the N'^LO and Mowfc potentials, respectively. 

In Fig. |4] bottom panel, we compare the components 
of the "exchange" -norm kernel up to ^ = 2. Contribu- 
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FIG. 5: (Color online. 



'Exchange" part of the diagonal norm 

1 -■ 
2 



kernel for the n-'*He(g.s.) | | channel as a function of the 



relative coordinates r and r', using the N^LO [U] (left) and 
Viowfc [28] (right) iViV potentials at hO. = 19 and 18 MeV, 
respectively. 





71 


72 


73 


N max 




Viowfe 




9 


-0.9547 


-0.06609 


-0.00310 


11 


-0.9539 


-0.06600 


-0.00288 


13 


-0.9530 


-0.06616 


-0.00290 


15 


-0.9526 


-0.06617 


-0.00292 


17 


-0.9524 


-0.06616 


-0.00293 


N 

^ ' max 




N^LO 




9 


-0.954 


-0.0633 


-0.00346 


11 


-0.945 


-0.0641 


-0.00452 


13 


-0.938 


-0.0643 


-0.00524 


15 


-0.933 


-0.0646 


-0.00599 


17 


-0.929 


-0.0645 


-0.00636 


19 


-0.927 


-0.0644 


-0.00661 


21 


-0.926 


-0.0645 


-0.00684 






AV14 




FY [30] 


-0.937 


-0.0663 


-0.00753 



TABLE I: The three largest negativ e ei genvalues of the "ex- 
change" part of the norm kernel (371 for the n-*He(g.s.) 
J^T =11 channel. Convergence with respect to the model- 
space size A'max of the NCSM/RGM results obtained using 
the Viowfe [28] and N^LO NN potentials at hQ = 18 and 19 
MeV, respectively. The calculated values for the AV14 NN 
potential of Ref. 00] are multiplied by -1 to adhere to the 
definition of the norm kernel adopted in the present paper. 
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FIG. 6: (Color on- 
line.) 

Diagonal n-*He(g.s.) 

potential 
kernels as a function 
of the relative coor- 
dinates r and r', us- 
ing the Mowfc |28] NN 
interaction. Model 
space and HO fre- 
quency are Nmax ~ 17 
and hO. = 18 MeV, re- 
spectively. 



eigenvalues are very close. The present results are also 
in good agreement (especially for N'^LO) with the eigen- 
values obtained in Ref. [SD] from a Faddeev-Yakubovsky 
calculation of the five-nucleon "exchange" norm, using 
the AV14 NN potential. 

The presence of a forbidden state affects also the po- 
tential kernels. The surface plots of Figs. [6]and[7]present 
"direct" and "exchange" potentials for the and 

I i channels, respectively. In the ^5*1/2 partial wave 
the Pauli-exclusion principle manifests itself again in the 
short-range repulsive action of the "exchange" potential, 
which effectively suppresses the interaction between one 
of the nucleons inside the a particle and the fifth nucleon, 
both in s-shell configuration. The situation is different 
in the channel, where the "exchange" kernel repre- 

sents a ~ 15% correction to the "direct" potential, and 
generates additional attraction. 

In the five-nucleon system the \^ \ ^s, the only forbid- 
den state (which is also the reason why the fivc-nuclcon 
g.s. occurs in P wave). For all other partial waves, the 
"exchange" part of the integral kernels introduces only 
a small deviation from orthogonality in the case of the 
norm, or small corrections to the "effective" n-a inter- 
action, in the case of the potential. These many-body 
corrections induced by the non-identical permutations in 
the inter-cluster anti-symmetrizers become less and less 
important with increasing relative angular momentum 
and have a limited range of about 5 fm. 



E. Orthogonalization 

The appearance of the norm kernel J^^iJ^{r',r) in 
Eq. ([4]) reflects the fact that the many-body wave func- 
tion ^''^ is expanded in terms of a non-orthogonal basis. 
Therefore, Eq. ^ does not represent a system of multi- 
channel Schrodinger equations, and "^(r) do not rep- 
resent Schrodinger wave functions. However, as we have 
seen in Sec. |IID 3| the non-orthogonality is short-ranged, 
as it originates from the non-identical permutations in 
the inter-cluster anti-symmetrizers. Thus, asymptoti- 



cally one has 



S{r' - r) 



(54) 



As a consequence the relative wave functions '^{r) 
obey the same asymptotic boundary conditions as the 
relative wave functions in a conventional multichannel 
collision theory, and it is possible to define physically im- 
portant quantities, such as, e.g., the scattering matrix, or 
the energy eigenvalues. The internal part of the relative 
wave functions, however, is still affected by the short- 
range non-orthogonality. Therefore, attention has to be 
paid when the latter wave functions are used to calculate 
further observables, such as, e.g., radiative capture cross 
sections, or, more in general, transition matrix elements. 

Alternatively one can introduce an orthogonalized ver- 
sion of Eq. (|4|, e.g., 

S{r' -r)- 
r'r 



E 



drr 



0. 



(55) 

where ]HI^,J^(r', r) is the Hermitian energy-independent 
non-local Hamiltonian defined by 



E 

7' 



dy'y' 



7 



dy 



X M,r^, {r\ y') 7i^,7(2/', y) U.J' {y, r) , (56) 

and the Schrodinger wave functions xl "^('') ^'^e the new 
unknowns of the problem, related to g^ "^('') through: 



E 



dyy^uhr.y)^-^^^ 

y 



(57) 



Here, Nl^,^{x',x) and N^,l{x' ,x) represent the square 
root and the inverse-square root of the norm kernel, re- 
spectively. In order to perform these two operations, we 
add and subtract from the norm kernel the identity in 
the HO model space 

r^^^^-' ^ r-^^-'-^) ^i?„,(,')i?„,(,) 



Rn't{r') A^J.^n Rni{r) 



(58) 
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FIG. 7: (Color on- 
line.) 

Diagonal n-'*He(g.s.) 
I" I (^-P3/2) potential 
kernels as a function 
of the relative coor- 
dinates r and r', us- 
ing the 1/iowfc 128] AfAf 
interaction. Model 
space and HO fre- 
quency are Nmax = 17 
and hQ. = 18 MeV, re- 
spectively. 



The matrix A"' ^ is the norm kernel within the truncated 
model space spanned by the HO Jacobi-channel states of 
Eq. (22 1. We give here the expression in the SNP basis 



[see also Eq. (37)]: 



A-' 



Pa- 



l.A 



(59) 

The generalization to the case of binary clusters with 
a > 1 is straightforward. Square root and inverse-square 
root of Af^,J^{r' , r) are then obtained by i) finding eigen- 
values, Ar, and eigenvectors, |(/9p of the matrix A"^ ^; 
a) calculating 



A ,^ 



E 



; (60) 



and, finally, in) replacing the model-space norm Aj^,,J 



in Eq. (58 1 with A^,„,^„ and A^,^, respectively, i.e., 



Y.Rn't{r')htl,unRni{r). 



(61) 



For the inverse operation to be permissible in Eq. ( 60 ) 



one has to exclude the subspace of (fully) Pauli-forbidden 
states for which Ar = (we note here that in the example 



of Sec. II D 3 the eigenvalues of the the norm kernel in 
the ^51/2 are related via Ar = 1 + 7r)- 

Both systems of coupled differential equations ^ 
and ( 55 1 can be cast in the form 



[T,.ei(r') + Vc(r')-(^- 



r It' u 



(r') 



J2 Jdrr^W;!,7ir',r) 



(r) 



= 0, 



(62) 



where 



(r) stands for either gi, (r) (in the non- 



orthogonalized case) or xi '^{f) (in the orthogonalized 
case), and W^,J'{r\ r) is the potential collecting all non- 
local terms present in the original equation. Obviously, 



in the (non-orthogonalized) case of Eq. ^ this non-local 
potential depends upon the energy. 

To provide some illustrative examples of non-local 
potentials corresponding to the orthogonalized case of 



Eq. ( 55 1 , we turn again to the n-a system, for which, as in 
Sec. | IID 3| we will present here results of single-channel 



calculations with the a particle in its g.s. Figure |8] shows 
the three partial waves ^S'1/2, and ^^3/2, obtained 

using the N^LO NN potential [^. The non local po- 
tentials for the three different spin-parity channels all 
rapidly vanish to zero beyond about 5 fm (as already ob- 
served in the non-orthogonalized integral kernels), while 
presenting substantially diverse structures at short range. 
We note in particular the strong repulsion between nu- 
cleon and a particle induced by the Pauli-exclusion prin- 



ciple in the 



1+1 



channel, and the potential well leading 



to the ^He resonance in the | 



I channel. 



F. Solution of the radial equation 



In solving Eq. ( 62 1 we assume that Vc (r) is the only in- 
teraction experienced by the clusters beyond a finite sep- 
aration To, thus dividing the configuration space into an 
internal and an external region. The radial wave function 
in the external region is approximated by its asymptotic 
form for large r. 



for scattering states, or 



(63) 
(64) 



for bound states. Here Hj (rj^, K^r) = Geir/^, K^r) =F 
^Fiiji^, K^r) are incoming and outgoing Coulomb func- 
tions, whereas Wi{r]^,K^r) are Whittaker functions. 
They depend on the channel state relative angular mo- 
mentum i, wave number d Sommerfeld parameter 
1]^. The corresponding velocity is denoted as v,^. The 
scattering matrix S^^ {i being the initial channel) in 



Eq. ( 63 1 , or binding energy and asymptotic normalization 
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8: (Color online. 



forthen-Q(g.s.) JT = | + i 



Orthogonalized non-local potentials 
I I and I I channels as func- 



tions of the relative coordinates r and r , usine the N^LO 
NN potential [27]. The index v - r " - - n+n. i i + i . 
replaced by the quantum number 



= {4g.s.0+0;l| + |;i^} is 
for simplicity. 



constant in Eq. ( 64 1 , together with the radial wave 



function in the internal region are obtained by applying 
to Eq. Q or to Eq. ( 55 1 the coupled-channel i?-matrix 
method on a Lagrange mesh [3T]. For the bound-state 
calculation depends on the studied binding energy. 
Therefore, the determination of the bound-state energy 
is achieved iteratively starting from an initial guess for 
the value of the logarithmic derivative of the wave func- 
tion at the matching radius rg. 

Finally, the accuracy of the i?-matrix method on a La- 



grange mesh is such that for a matching radius of tq = 15 
fm, N = 25 mesh points are usually enough to determine 
a phase shift within the sixth significant digit. The typi- 
cal matching radius and number of mesh points adopted 
for the present calculations are rp = 18 fm and N = 40. 



III. RESULTS 
A. A = 4 

The four-nucleon scattering problem, with its compli- 
cated interplay of low-energy thresholds and resonances, 
represents a serious theoretical challenge, only recently 
addressed by means of accurate ab initio calculations. 
Important developments in the numerical solution of the 
four-nucleon scattering equations in momentum space [5] , 
and in the treatment of the long-range Coulomb interac- 
tion [32| have led to very accurate ab initio calculations 
of scattering observables in the energy region below the 
three-body break-up threshold. 

In this section we use the four-nucleon system as a test- 
ground to study the performances of our newly-developed 
NCSM/RGM approach within the SNP basis. In partic- 
ular, we present here results of coupled-channel calcula- 
tions restricted to basis channel states with the three- 
nucleon target in its g.s. (corresponding to channel in- 
dexes of the type v = {3 g.s. | 5; ^] si}). Indeed, 
we are interested to the energy region below the break-up 
threshold of the A = 3 target. 

We start by studying the convergence of our calcula- 
tions with respect to the HO model-space size (A^max) for 
the simplest of the A = 4 scattering channels, i.e., the n- 
•^H. This is a purely T — 1 system, with no Coulomb 
interaction between target and projectile. As the over- 
all convergence behavior strongly depends on the model 
of NN interaction adopted, we first consider results ob- 
tained using the "bare" Viowk potential [2E1- These are 



summar ized in Table |ll| Both g.s. energy and n-^H 
scattering data present a rather weak dependence on 
Nmax- However, a sudden worsening in convergence rate 
is noticeable in the higher model spaces, especially for the 
phase shifts of small magnitude. This is in part a reflec- 
tion of the sharp cutoff-function used to derive the Viowfe 
potential (here we use the version derived from AV18 
with cutoff A = 2.1 fm-i). 

Next we present n-^H phase shifts obtained using the 
N'^LO NN interaction [27]. The convergence behavior 
shown in Fig. [9] was achieved using two-body effective 
interactions tailored to the model-space truncation, as 
outlined in Sec. |IIB For the ^Sq, ^Pi and ■^S'l partial 
waves, the increase in model-space size produces gradu- 
ally smaller deviations with a clear convergence towards 
the A^max — 19 results. The rest of the phase shifts, par- 
ticularly the ^Pq, show a more irregular pattern. Nev- 
ertheless, in the whole energy-range we find less than 2 
deg absolute difference between the phases obtained in 
the largest and next-to-largest model spaces. The agree- 
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TABLE 11: Calculated g.s. energy (in MeV) and n-^H phase shifts (in degrees) and total cross section (in barns) for 
increasing A^max at HQ = 18 MeV, obtained using the Viow* AA'^ potential [2H]- The scattering results were obtained in a 
coupled-channel calculation including only the g.s. of the nucleus (i.e. the channels v = {3g.s. ^ \\ l|^|i s^}). 











n-^'H (Skin = 0.40 MeV) 






^ max 


-Eg.s. 


0+ (^So) 


0- (^Po) 


1+ QSx) 1- (iPi) 1- (^Pi) 


1- (.) 2- (^Pa) 





9 


-7.80 


-20.2 


0.93 


-18.9 


0.85 


1.96 


-18.0 


3.01 


0.99 


11 


-7.96 
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1.15 


13 
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15 
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1.31 


17 


-8.12 


-25.1 


1.06 


-22.6 


1.52 
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n-^H (E'kin = 0.75 MeV) 
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FIG. 9: (Color online) Calculated n-"^H phase shifts as a function of the relative kinetic energy in the cm. frame i?kin, using the 
N^LO A'A^ potential [27| in the model spaces A^max = 11 — 19, at HQ, = 22 MeV. All results were obtained in a coupled-channel 
calculation including only the g.s. of the "^H nucleus (i.e. the channels v — {3g.s. | |; l| |; s£}). 



ment within 1.5 deg of the A^max = 19 results obtained 
with two different HO frequencies, Ml = 19 and hil — 22 



MeV, (see Fig. 10) is a further indication of the fairly 



good degree of convergence of our calculation. 

In order to verify our approach, in Fig. 10 we com- 



pare our n-^H results to earlier ab initio calculations 
performed in the framework of the Alt, Grassberger and 
Sandhas (AGS) equations [H1IS3], using the same N'^LO 
NN potential. We note that in general the agreement 
between the two calculations worsens as the relative ki- 
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FIG. 10: (Color online) Calculated 
n-^H phase shifts using the N"^LO 
NN potential [57] for iV^ax = 19 
and hn = 19, and 22 MeV, com- 
pared to AGS results of Refs. pi33|. 
All NCSM/RGM results were ob- 
tained in a coupled-channel calcula- 
tion including only the g.s. of the 
■^H nucleus (i.e. the channels u = 
{3g.s.i + i;li+|;.n). 



FIG. 11: (Color onhne) Calculated 
n-^H phase shifts using the CD- 
Bonn iViV potential [34] for iVmax = 
19 and HQ — 19 MeV, compared 
to AGS results of Refs. [H [33]. 
All NCSM/RGM results were ob- 
tained in a coupled-channel calcula- 
tion including only the g.s. of the 
■^H nucleus (i.e. the channels v = 
{3g.s.| + i; ir|;.^}). 



netic energy in the cm. frame, -Ekin, increases. For the 
P-waves in particular we can reasonably reproduce the 
AGS calculation for energies within 1 MeV while we can 
find differences as large as 17 deg ('^P2) at ii^kin = 2.6 
MeV. In Fig. [11] an analo gous comparison performed for 
a second realistic NN interaction, the CD-Bonn poten- 
tial [34], leads to a similar picture. (Note that, as for 
N^LO, the NCSM/RGM results for CD-Bonn were also 
obtained using two-body effective interactions.) These 
discrepancies are due to the influence, increasing with 
energy, played by closed channels not included in our 
calculations, such as those with the A — 1 — 3 eigenstates 
above the I^^ = g.s., and {A — a~2, a — 2) configu- 
rations, present in the AGS results. As an indication, in 
Ref. 2] it was shown that the omission of three-nucleon 



partial waves with ^ < /i < | leads to effects of compa- 
rable magnitude on the AGS results, especially for the 
351,3 Pi and 3P2. 

All A— I = 3 states but the I^^ = g.s. are in the 
continuum, and correspond to a break-up of the three- 
nucleon target. Therefore, the corresponding {A — a — 3, 
a—1) channels do not represent "open" rearrangement 
channels in the energy range considered here. However, 
it is clear from the previous analysis that the virtual ex- 
citation of the A— 1 = 3 target has an important influence 
on the n-^H elastic phase shifts, and should be included 
in the NCSM/RGM approach in order to reach full con- 
vergence, and hence agreement with the AGS calcula- 
tion. Obviously, considering the localized nature of the 
NCSM wave functions, for each /J^^ 7^ 5"*" one obtains a 
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large series of positive-energy eigenstates corresponding 
to a denser and denser discretization of the A — l = 3 
continuum, as the HO models space increases. Conse- 
quently, it would not be conceptually sound to try and 
include these states in the NCSM/RGM SNP basis, not 
to mention that it would not be computationally fea- 
sible either. On the other hand, the A = 4 low-lying 
spectrum contains a finite number of fairly narrow reso- 
nances, which can be reasonably reproduced diagonal- 
izing the four-body Hamiltonian in the NCSM model 
space. Therefore, it is clear that the most efficient way of 
tackling the A — 4 scattering problem would be for us to 
use an over-complete model space formed by both tradi- 
tional NCSM four-body states and NCSM/RGM cluster 
states. Although it is in our intentions to pursue this 
approach, we leave it for future investigation. 

In the remaining part of this Section we will discuss 
the scattering of protons on '^He targets. This is once 
again a purely T = 1 system, but differs from the n-^H 
case because of the presence of the Coulomb interaction 
between the clusters, both charged. The treatment of 
the Coulomb interaction between target and projectile, 
as explained in Sec.|ll] does not represent a major obsta- 
cle in the NCSM/RGM approach. In particular, in the 
following we will show that the p-^He phase shifts present 
a similar convergence trend as the one observed in their 
neutral counterparts. 



In order to perf orm a direct comparison with the n- 
data, in Table 



III 



we present He g.s. energy and 
p-'^He scattering phase shifts for the same ("bare") Mowfe 
NN potential [2S] and relative kinetic-energy values as 
in Table [III expected, the growth of the nuclear phase 
shifts from the zero energy is slower in the presence of 
the Coulomb repulsion between the clusters. This is es- 
pecially visible at the very low energies considered here 
(-E-kin = 0.4, and 0.75 MeV). As the scattering data, par- 
ticularly in the P waves, are very small in magnitude, the 
somewhat slower convergence rate in the biggest model 
spaces already noticed in the n-^H case is emphasized 
even more here. This feature, partly related to the sharp 
cutoff of the Viowk potential, results in differences of a few 
tenths of a degree between the A^max — 17 and A^max = 19 
phase shifts. 



Figure [l2| shows p-'^Hc phase shifts obtained using the 
N'^LO A^A^ potential, the results of this work (solid and 
dashed lines) and those of AGS calculations [33^ (-I-). 
The use of two-body effective interactions tailored to the 
size of the adopted model-spaces, guarantees also in this 
case a fairly good agreement (of the same order as in 
Fig. [9]) between the A^ax — 17 and A^ax = 19 calcula- 
tions. The comparison to the AGS results shows that the 
NCSM/RGM SNP basis with the ^He nucleus in its g.s. 
provides the bulk of the p-^He elastic phase shifts, con- 
firming the observations made for the n-'^H scattering. 



B. 



Driven by wider efforts to develop a predictive ab initio 
theory of low-energy reactions on light nuclei, ab initio 
calculations for scattering processes involving five nucle- 
ons are beginning to be realized in the last couple of 
years, but are still a rare exception. First, the n-a low- 
lying = 3/2~ and l/2~ P-wave resonances as well as 
the 1/2+ 5- wave non-resonant scattering below 5 MeV 
cm. energy were obtained using the AV18 A^A^ potential 
with and without the three- nucleon force, chosen to be ei- 
ther the Urbana IX or the Illinois-2 model . The results 
of these Green's function Monte Carlo (GFMC) calcula- 
tions revealed sensitivity to the inter-nucleon interaction, 
and in particular to the strength of the spin-orbit force. 
Soon after, the development of the ab initio NCSM/RGM 
approach allowed us to calculate both n- and (for the first 
time) p-a scattering phase shifts for energies up to the 
inelastic threshold, using realistic A A potentials [T7]. In- 
deed, nucleon-a scattering provides one of the best-case 
scenario for the application of the NCSM/RGM approach 
within the SNP basis. This process is characterized by 
a single open channel up to the *He break-up thresh- 
old, which is fairly high in energy. In addition, the low- 
lying resonances of the *He nucleus are narrow enough 
to be reasonably reproduced diagonalizing the four-body 
Hamiltonian in the NCSM model space. Therefore, they 
can be consistently included as closed channels in the 
NCSM/RGM SNP model space. In the following we give 
a detailed description of previously published |17j and 
new results for nucleon-a scattering. 

First we present single channel calculations carried 
out using n-a channel states with the a particle in 
its g.s., i.e., characterized by the channel index ly = 
{4 g.s. 0+0; 1 1"*"!; I (or simply by the angular quan- 



tum number £). In particular. Table IV shows the good 
degree of convergence with respect to A^ax obtained for 
the *He g.s., and for the n-a (^S'1/2, and ^^3/2) 

phase shifts and total cross section at i?kin = 2.5 and 5 
MeV, using the (bare) Mowfe A^A^ interaction. The corre- 
sponding p-a scattering phase shifts can be found in Ta- 
ble |V] The HO model-space dependence of the Viowfc n-a 
phase shifts is presented also in the left panel of Fig. [13] 
where it is explored for a wider range of energies, and 
compared to an analogous plot for the N'^LO AA^ inter- 
action (central panel) . Despite the use of two-body effec- 
tive interaction as outlined in Sec. [II B[ the convergence 
rate is visibly much slower for N^^LO. This gives a mea- 
sure of the stronger short-range correlations generated by 
this potential. The ^P3/2 phase shifts present the largest 
(up to 5 deg in the energy range between 1 and 4 MeV) 
differences between the N^nax = 15 and 17 calculations, 
which are otherwise no more than 2 deg apart. 



The third (right) panel of Fig. 13 compares the A'„iax = 
17 results for the previously discussed Mowfc and N'^LO 
A^A^ interactions, and those obtained with the CD-Bonn 
NN potential [SJ. The NCSM/RGM calculations for the 
latter potential were carried out using two-body effective 
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TABLE III: Calculated ^ 


He g.s. energy 


(in MeV) and p-' 


*He phase shifts (in degrees) for increasing 


A"max at HQ = 


18 MeV, 


obtained 


using the V[owfc 


A^A'' potential 


|28|. The scattering results were obtained in a coupled-channel calculation 


including 


only the 


g.s. of the ^He nucleus (i.e. the channels u — {3g 


s i + i 

■°- 2 2 


l| + |;sn). 
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1 c: 
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FIG. 12: (Color onhne) Calculated 
p-^He phase shifts for the N^LO NN 
potential in the model spaces 
Af^ax = 17 - 19, at no. = 22 MeV, 
compared to AGS results of Ref. |33| . 
All NCSM/RGM resufts were ob- 
tained in a coupled-channel calcula- 
tion including only the g.s. of the 
■^He nucleus (i.e. the channels v = 
{3g.s.i + i;li+|;.n). 



interactions, and present a convergence pattern similar 
to the one observed for N'^LO. Clearly, the and 
phase shifts are sensitive to the interaction models, 
and, in particular, to the strength of the spin-orbit force. 
This observation is in agreement with what was found 
in the earlier study of Ref. [9 . Following a behavior 
already observed in the structure of p-shell nuclei, CD- 
Bonn and N'^LO interactions yield about the same spin- 
orbit splitting. On the contrary, the larger separation 



between the Mowfc 3/2~ and l/2~ resonant phase shifts 
is direct evidence for a stronger spin-orbit interaction. 

As the 1/2+ channel is dominated by the repulsion 
between the neutron and the a particle induced by the 
Pauli exclusion principle (see also Sec. IID3I, the short- 
range details of the nuclear interaction play a minor role 
on the ^5*1/2 phase shifts. As a consequence, we find very 
similar results for all of the three adopted NN potential 
models. Worth of note is the different behavior of the 
CD-Bonn results close to the zero energy, which appears 
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FIG. 13: (Color online) Dependence on Nmi^x of the n-a{g.s.) phase shifts with the Viowfc i28j (left panel) and N'^LO ^\ (central 
panel) A'^A'^ potentials at hil = 18 and 19 MeV, respectively. In the right panel, results obtained in the largest model space 
(A^max = 17). The calculation for the CD-Bonn ^ NN interaction was performed at hD, = 19 MeV. 



TABLE IV: Calculated *He g.s. energy (in MeV) and n-*He 
phase shifts (in degrees) and total cross sections (in barns) for 
increasing A^max at HQ = 18 MeV, obtained using the Viowfc 
A'^A'^ potential |28) . The scattering results were obtained in a 
single-channel calculation including only the g.s. of the *He 
nucleus (i.e. the channel u — {4g.s.0"'"0; l|^|; |^})- 
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17 






-58.6 
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also in the P waves. 

Next we explore the effect of the inclusion of excited 
states of the *He on the n-a scattering phase shifts ob- 
tained with the N'^LO NN interaction. In contrast to 
the A = A scattering, discussed in the previous section, 
binary channels of the type {A— 2, 2) have here a much 
suppressed effect due to the large binding energy of the 
'^He nucleus. However, in order to reach full convergence 
it is still necessary to take into account the virtual exci- 
tations of the A—\ — 4 target. To this aim we extend the 
NCSM/RGM SNP model space to include closed chan- 



nels of the type v ^ {A P'ex. Il^Ti, \; si} with 
/^^iTi = 0+0, 0^0, 1-1, 2-0, and 2-1, and "P*ex." spec- 
ifies that, for each of these spin-parity and isospin com- 
binations, we consider only the first (low- lying) excited 
state. 

In addition to the above discussed single-channel re- 
sults (dotted line). Figure 14 shows coupled-channel cal- 
culations for five different combinations of ''He states, 
i.e., i) g.s. ,0+0 (dash-dotted line), ii) g.s. ,0+0, 0"0 (dash- 
dot-dotted line). Hi) g.s.,0+0, 0"0, 1"0, 1"I (dash-dash- 
dotted line), iv) g.s.,0+0, 2-0 (dashed line), and v) 



TABLE V: Calculated p-^He phase shifts (in degrees) for in- 
creasing A'^max at hQ, — 18 MeV, using the Mowfc NN poten- 
tial [2^. Results were obtained in a single-channel calculation 
including only the g.s. of the ^He nucleus (i.e. the channel 
;/ = {4g.s.0+0; li+|; ^£}). 
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FIG. 14: (Color online) Influence of the lowest six excited states {I^^Ti = 0+0, 0~0, 1~0, 1~1,2"0,2~1) of the a particle on 
the n-a (left panel), (central panel), and ^^3/2 (right panel) phase-shift results for the N'^LO NN potential [27] 

at hfl — 19 MeV. Dotted (g.s.) and dash-dotted (g.s., 0'''0) lines correspond to single- and coupled-channel calculations in a 
Nmax ~ 17 model space, respectively. The effects on the and phase shifts of the further inclusion of, respectively, 

the 0~0, 1~0, 1~1, and 2~0, 2~1 states are investigated in a N^^,^ — 15 model space. 



g.s., 0+0, 2^0, 2^1 (solid line). The 0+0 excited state has 
a minimal influence on all three phase shifts. In addition, 
for ^5*1/2 (left panel) no further corrections are found 
in the four larger Hilbert spaces obtained by including 
the low-lying negative parity states of ^He (for clarity of 
the figure we omitted these latter results). On the 

contrary, we find larger deviations on the (central 
panel) and ^^3/2 (right panel) phase shifts, after inclu- 
sion of the 0^0, 1~0, and 1~1 states for the first, and 
of the 2~0 and 2~1 states for the second. These nega- 
tive parity states influence the P phase shifts, because 
they introduce couplings to the s-wave of relative mo- 
tion. Though also I^^ = 1~ couples to £ = in the 3/2^ 
channel, the coupling of the I^^ — 2~ states is dominant 
for the ^^3/2 phase shifts. 

Figure [15] provides further evidence that the 
NCSM/RGM SNP model space formed by nucleon-a bi- 
nary channels with the a particle in its ground and first 
0+0 excited states is sufficient to reach full convergence of 
the ^51/2 phase shifts, also in presence of the Coulomb re- 
pulsion between proton and a particle. In the left panel, 
both n- and p-a N'^LO results show negligible depen- 
dence on the HO frequency, when varied from hfl — 19 
to 22 MeV. In the right panel, the latter phase shifts and 
the corresponding Viowfc and CD-Bonn ^S'1/2 results are 
compared to an accurate multichannel i?-matrix analy- 
sis of nucleon-a scattering. The overall best agreement 
with experiment (quite remarkable for p-a) is obtained 
for the CD-Bonn NN interaction, where the different be- 
havior of this potential near the zero energy is favored by 
the data. The N'^LO phase shifts are not very dissimilar, 
and reproduce the i?-matrix analysis starting from an en- 
ergy of roughly 2 MeV. The Viowfe interaction generates 
the largest deviation from experiment. While these are 
"residual" reflections of the interaction details, otherwise 
masked by the Pauli exclusion principle, it becomes ev- 
ident that scattering calculations can provide important 
additional constrains on the nuclear force. 



A comparison to the i?-matrix analysis of Ref. |35j 



including Pi 
sented in Fig. 16 



2 p 

/2j -'3/21 



and D^i2 partial waves, is pre- 



Here, the n- (left panel) and p-a 
(right panel) phase shifts were obtained with the N'^LO 
NN potential, including the first six ^He excited states, 
as shown in Fig. 14 The magnitude of the ^Z?3 



.-3/2 phase 

shifts, calculatedTas the in a NCSM/RGM SNP 

model space with ground and first 0+0 excited states of 
the a particle, is qualitatively reproduced. On the con- 
trary, the P phase shifts present both insufficient magni- 
tude and splitting with respect to the predictions of the 
i?-matrix analysis. Although the inclusion of two more 
■^He negative excited states (first 0^1 and second 1~1) 
beyond the five considered here could introduce small 
corrections, it is not likely that they would explain the 
present discrepancy with respect to experiment of the 



'P^ 



1/2 



and ^^3/2 results. 



On the other hand, consider- 



ing the sensitivity of these phase shifts to the strength 
of the spin-orbit force, the inclusion of the NNN terms 
of the chiral interaction would probably lead to an en- 
hanced spin-orbit splitting, and recover the predictions 
of the i?-matrix analysis. 
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With the advent of experimental programs on exotic 
nuclei, the description of weakly bound nuclei has be- 
come one of the priorities of modern nuclear theory. As 
techniques traditionally successful for well-bound nuclei 
struggle to reproduce new phenomena observed in the 
radioactive-beam facilities, the interplay of structure and 
reaction mechanisms is now unanimously recognized as 
a prime element for a successful description of weakly 
bound nuclei. Such interplay is an intrinsic characteristic 
in the ah initio NCSM/RGM, where bound and scatter- 
ing states are treated in a unified formalism. In this Sec- 
tion we test the performance of our formalism in the SNP 
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FIG. 15: (Color online) A'^-q ^51/2 phase shifts as a function of the relative kinetic energy in the cm. frame -Bkin- Frequency 
dependence (left panel), and comparison to an i?-matrix analysis of data (+) 35 (right panel) of coupled-channel calculations 
including both ground and first 0+0 excited states of the a particle, in a Nmax = 17 model space. 
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FIG. 16: (Color online) Calculated n- (left panel) and p-a (right pa nel) phase shifts for the N^LO NN potential [27] compared 
to an i?-matrix analysis of data (-I-) [3 
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The ^Pi 



12 and ^P3/2 phase shifts correspond to the 



Si/ 2 results as in Fig 
dash-dash-dotted (g.s.,0+0, 0"0, 1"0, 1"1) and solid (g.s.,0+0, 2^0, 2~) lines of Fig 

in a coupled-channel calculation including ground and first 0+ excited state of ''He, in a Nmux = 17 HO model space 
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The -Da/2 phase shifts were obtained 



basis for the description of one-nucleon halo systems, and 
at the same time show the versatility and promise of the 
NCSM/RGM for the description of the structure and re- 
actions of p-shell nuclei. 

Among light drip-line nuclei, ^^Be provides a conve- 
nient test of several important properties of neutron rich 
nuclei. In particular, the parity-inverted ground state of 
this nucleus, first observed by Talmi and Unna in the 
early 1960 's represents one of the best examples of 
disappearance of the N = 8 magic number with increasing 
N/Z ratio. 

The only previous ab initio investigations of the ^^Be 
low-lying states, consisting of large-scale NCSM calcula- 
tions with realistic NN potentials, were unable to repro- 



duce this phenomenon f37^. This result was partly at- 
tributed to the size of the HO basis, which was not large 
enough to reproduce the correct asymptotic of the n-^°Be 
component of the 11-body wave function. At the same 
time the calculations performed with the INOY (inside 
non-local outside Yukawa) NN potential of Doleschall et 
al. [38] suggested that the use of a realistic NNN force in 
a large NCSM basis might correct this discrepancy with 
experiment. 

The correct asymptotic behavior of the n-^^Be wave 
functions can be reproduced when working within micro- 
scopic cluster techniques. Starting from a microscopic 
Hamiltonian containing the Volkov NN potential [55] . 
the Coulomb interaction, and a zero-range spin-orbit 
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T — I — I — I — I — I — I — I — I — I — r 

1° Be 



'Be states 

g.s.,2t,2j,l| 

-- g.s.,2+2j 

. . . g.s. 




mg on ^°Be and low-lying levels of ^^Be are stud- 
ied by means of NCSM/RGM coupled channel calcula- 
tions with n-^°Be channel states including ^°Be ground, 
and 2^,2^,1]'' excited states (corresponding to chan- 
nel indexes of the type v = {Wa\Tl^T\\ l|^|; s^} 
with axIl^Tx= g.s. 0+1, P*ex.2+1, 2"'^ex.2+l, and 
P*ex. 1+1). The NCSM i°Be eigenstates, calculated for 
the first time in occasion of the publication of Ref . ^42] , 
are obtained here in a A^max = 6 model space. Cor- 
respondingly, the 11-body NCSM/RGM model space is 
.^max = 6(7) for negative-(positive-)parity wave functions. 
In order to perform a direct comparison to the NCSM re- 
sults for -'^^Be [37] obtained using the CD-Bonn NN in- 
teraction |34j . we adopt the same nuclear potential and 
optimal HO frequency, h£L = 13 MeV. 

we present energies of the lowest 1/2+ 
states of ^^Be obtained in the NCSM and in 
the NCSM/RGM calculations. Clearly, there is httle dif- 
ference between the TV^ax = 6/7 and A'^max = 8/9 NCSM 
results. The l/2~ state is the ground state and the exci- 
tation energy of the 1/2+ state is about 3 MeV (or about 
2.8 MeV above the n-i°Be threshold). The NCSM/RGM 
calculations that include ^"Be g.s. and the three lowest 
calculated excited states (2]^, 2^, and 1+) show little 
change for the energy of the 1/2" state. However, we 
observe a dramatic decrease (~3.5 MeV) of the energy of 
the 1/2+ state. In the NCSM/RGM calculations, both 
the l/2~ and 1/2+ states are bound and the 1/2+ state 
becomes the ground state of ^^Be. To understand the 
binding mechanism of the 1/2+ state, we evaluated mean 
values of the relative kinetic and potential energies as well 
as the mean value of the ^"Be energy. The resuls are given 
in Table |VII| together with the corresponding values ob- 
tained by restricting all the kernels to the model space. 
The modcl-space-restricted calculation is then similar, al- 
though not identical, to the standard NCSM calculation. 
In particular, we lose the correct asymptotic behavior of 
the wave functions guaranteed in the NCSM/RGM. We 
observe that both the relative kinetic energy and poten- 
tial energies decrease in the full NCSM/RGM calculation. 
The drop of the relative kinetic energy is significantly 
more substantial than of the potential energy, due to the 
re-scaling of the relative wave function in the internal re- 
gion when the Whittaker tail is recovered. This is the 
main cause of the dramatic decrease of the energy of the 
1/2+ state, which makes it bound and even leads to a 
g.s. parity inversion. Although we cannot exclude that 
the NNN force plays a role in the inversion mechanism, 
it is clear that a proper treatment of the coupling to the 
ri-^'^Bc continuum is essential in explaining the g.s. par- 
ity inversion. 

Our calculated ■^S'1/2 n-^°Be phase shifts are displayed 
in Fig. [17] We show results obtained with different num- 
ber of ^'Be states. The phase shift does not change sig- 
nificantly, once the lowest 2+ state is taken into account. 
A bound state was found, however, already by using just 
the ^°Be ground state. We also calculated the S-wa,ve 
scattering length. With all four ^"Be states (g.s., 2^ , 2^, 



1 1.5 2 

iJkin [MeV] 

FIG. 17: (Color online.) Calculated ^5*1/2 n-^°Be phase 
shifts as a function of Skin, using the CD-Bonn A'' A*' poten- 
tial. NCSM/RGM calculation as in Table [yi] The obtained 
scattering length is -1-10.7 fm. 



force [40) . Descouvemont was able to reproduce the in- 
version of the 1/2+ and 1/2" ^^Be bound states within 
the generator coordinate method (GCM) [U]. However, 
the use of two different parameterizations of the Volkov 
potential for positive- and negative-parity states (chosen 
to reproduce, respectively, the experimental binding en- 
ergies of the 1/2+ g.s., and 1/2" first excited state) was 
key to this result. With a single parametrization for both 
parities, the lowest energy is obtained once again for the 
1/2" state, in contradiction with experiment. The in- 
troduction of the tensor force (missing in ^5T]) and the 
use of a richer structure for the ^^Be wave function could 
probably cure this problem. 

A more complete bibliography on the ^^Be g.s. parity- 
inversion and the theoretical attempts to reproduce it 
can be found in Refs. [37], [H], and references therein. 

Here, low-energy phase shifts for neutron scatter- 
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FIG. 18: (Color online.) Calculated ^Ds/a n-^"Be phase 
shifts as a function of iSkin, using the CD-Bonn A^TV potential. 
NCSM/RGM calculation as in Table jVl] 
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TABLE VI: Calculated energies (in MeV) of the ^"Be g.s. and of the lowest negative- and positive-parity states in ^^Be, obtained 
using the CD-Bonn NN potential 01] at ft^ = 13 MeV. The NCSM/RGM results were obtained using n-|-^°Be configurations 
with A'^max = 6 g.s., 2^, 2^, and 1^ states of ^"Be. 







i«Be 


"Be(r 
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"Be(r 
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max 


Eg.s. 


E 


Eth 


E 


Eth 


NCSM [3Zl|42] 


8/9 


-57.06 


-56.95 


0.11 


-54.26 


2.80 


NCSM [13111]," 


6/7 


-57.17 


-57.51 


-0.34 


-54.39 


2.78 


NCSM/RGM" 






-57.59 


-0.42 


-57.85 


-0.68 


Expt. 




-64.98 


-65.16 


-0.18 


-65.48 


-0.50 



"present calculation 



and 1^) included we found the ^5*1/2 scattering length of 
-1-10.7 fm. This can be compared to the value of -1-13.6 
fm obtained in the GCM calculations of Ref. [H]. In 
those calculations, the ^^Be experimental binding energy 
of 0.5 MeV was fitted. Our calculated binding energy is 
slightly higher: 0.68 MeV. Correspondingly, our calcu- 
lated scattering length is smaller. 

In Fig. 18 we show our calculated ^D^j2 n-^°Be phase 
shifts obtained using different number of ^"Be states. We 
find a resonance in this channel below 3 MeV. To observe 
this resonance, it is crucial to include at least the first 2+ 
excited state of ^^Be in the NCSM/RGM calculations. 
A restriction to just the ^°Be ground state results in a 
smooth slowly rising phase shift with no resonance at 
low energy. We note that in experiment, there is a res- 
onance at ~ 1.8 MeV with a tentative spin assignment 
(5/2,3/2)+ m- 

Results presented in this subsection demonstrate the 
promise of the NCSM/RGM approach for applications 
to the p-shell nuclei. A significant improvement in the 
description of halo nuclei is achieved in particular when 
comparing to the standard NCSM calculations. The "Be 



and 



'Be calculations discussed here were obtained in 



a limited model space. We will improve on this in the 
future by expanding the model space sizes similarly as 
we did for the A = 4 and A = 5 systems. This will allow 
us to calculate reliably also e.g. the P-wave scattering 
length predicted to have a very large value [H]. For light 
p-shell nuclei, it is presently feasible to perform calcula- 
tions with A^max ~ 12 — 16. For heavy and mid-p-shell 
nuclei, it becomes possible to employ the importance- 
trunctated NCSM [IS] to expand target wave functions 
in large A'max model spaces. 



TABLE VIL Mean values of the relative kinetic and potential 
energy and of the intenal ^"Be energy in the ^^Be 1/2+ ground 
state. AU energies in MeV. NCSM/RGM calculation as in 
Table IVll See the text for further details. 



NCSM/RGM 




(W) 


iJrBe(g.s.,ex.)] 


-Etot 


Model Space 


16.65 


-15.02 


-56.66 


-55.03 


Full 


6.56 


-7.39 


-57.02 


-57.85 



IV. CONCLUSIONS 

We have presented in detail the NCSM/RGM formal- 
ism. This is a new ab initio many-body approach capa- 
ble of describing simultaneously both bound and scatter- 
ing states in light nuclei, by combining the RGM with 
the use of realistic interactions, and a microscopic and 
consistent description of the nucleon clusters, achieved 
via the ab initio NCSM. In particular, we have derived 
the algebraic expressions for the integral kernels for the 
case of a singlc-nucleon projectile, working both with 
the Jacobi-coordinate, and SD single-particle coordinate 
bases. As the spurious cm. components present in the 
SD basis were removed exactly, in both frameworks the 
calculated integral kernels are translationally invariant, 
and lead to identical results. Several analytical as well 
as numerical tests were performed in order to verify 
the approach, particularly by benchmarking independent 
Jacobi-coordinate and SD calculations for systems with 
up to 5 nucleons. 

Among the applications, we presented results for neu- 
tron scattering on ■^H, ''He and ^"Be and proton scatter- 
ing on '^'^He, using realistic A^A^ potentials. Our A — A 
scattering results were compared to earlier ab initio cal- 
culations. We found that the CD-Bonn NN potential in 
particular provides an excellent description of nucleon- 
"'He S'-wave phase shifts. On the contrary, the P-wave 
phase shifts that we obtained with any of the realistic 
NN potentials present both insufficient magnitude and 
splitting with respect to the P-matrix analysis of the 
data. It is anticipated that the inclusion of the NNN 
terms of the chiral interaction would lead to an enhanced 
spin-orbit splitting, and recover the predictions of the R- 
matrix analysis. An important topic of this work has 
been the investigation of the parity inversion of the -'^Be 
nucleus. Although we cannot exclude that, e.g. the 
NNN force plays a role in the inversion mechanism, 
we have demonstrated that a proper treatment of the 
coupling to the n-^'^Be continuum leads to a dramatic 
decrease of the energy of the state, which makes it 
bound and even leads to a g.s. parity inversion. 

It is straightforward to extend the NCSM/RGM for- 
malism to include two- nucleon (deuteron) , three-nucleon 
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(triton and ^He) and four-nucleon (^He) projectiles. Fur- 
ther, it is possible and desirable to extend the binary- 
cluster (A-a, a) NCSM/RGM basis by the A-nucleon 
NCSM basis to unify the original ab initio NCSM and 
the NCSM/RGM. In this way, a much faster convergence 
of many-body calculations will be achieved compared to 
the original approaches and, most importantly, an opti- 
mal and balanced unified description of both bound and 
unbound states will be obtained. In the NCSM/RGM 
a large HO basis expansion is needed not just for con- 
vergence of the target and projectile eigenstates but also 
for convergence of the localized parts of the integration 
kernels. The recently developed importance-truncated 
NCSM [35 makes use of large A/'max model spaces possi- 
ble even for heavy p-shell nuclei and beyond. A utility of 
the importace-truncated target wave functions within the 
NCSM/RGM formalism will allow convergence of scat- 
tering calculations for heavier nuclei similarly as it was 
demonstrated here for A — A and A — b systems. Devel- 
opment of the two-to-four-nucleon projectile formalism, 
the unification of the NCSM/RGM with the standard 
NCSM, that we name the ah initio NCSM with con- 
tinuum (NCSMC), and applications of the importance- 
truncated wave functions within the NCSM/RGM are 
under way. 



Acknowledgments 



LDRD grant PLS-09-ERD-020, and from the U. S. 
Department of Energy Grant DE-FC02-07ER41457 is 
acknowledged. 



APPENDIX A: JACOBI-COORDINATE 
DERIVATION 



Continuing from Sec. II D 1 where we have discussed 
the exchange part of the norm kernel, here we complete 
the Jacobi-coordinate derivation of the integral kernels 
for the A — ?> [a—l) system. For the notation we refer 



(461 



the interested reader to Eqs. (44l- 

As shown in Eq. (40 1, in the case of the "direct" - 



potential kernel one needs to evaluate matrix elements 
of the interaction between the last two nucleons, {V{r2 — 
cr2cr3r2r3)(l — -P23)). It is therefore convenient to in- 
troduce two new Jacobi coordinates 



Ci 

C2 



2 rl 



{r2 + ^3) - ri 



V2 



^2 -ra 



(Al) 
(A2) 
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and switch to the HO basis states in which nucleons 2 
and 3 are coupled together to form two-particle states 
of the form (C20'2cr3T2T3|A^2i2'S'2>/2T2), where N2,L2 are 
the HO quantum numbers corresponding to the harmonic 
oscillator associated with ^2, and S'2, J2, and T2 are the 
two-nucleon spin, total angular momentum and isospin 
quantum numbers of the (2,3)-nucleons couple, respec- 
tively. This task can be achieved, e.g., continuing from 
the expansion of Eq. ( 46 1 : 



^20-1 0-2 0-3 



(ni^i,n£)A; {s,-)Z 



2 2 ''1 



5 Z S: 



{N2L2,MC,A\ni,ni£i,A). 

2 ) N2L2MC 



2 - ^2 ) S2 
C \ J 

xJ2LZJJ2{ L2 S2 J2 \ (ClC2<Jia2a3nT2T3 [NCJ-N2L2S2J2T2]J^T) . 

A Z J 



(A3) 



Here A/", C, and are the HO quantum numbers corre- tal angular momentum of the first nucleon with respect to 
sponding to the HO state associated with Ci, and the to- the center of mass of the last two, respectively. Further, 
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{N2L2,AfC,A\ni,ni£i,A)3 are the general HO brackets 
for two particles with mass ratio 3, which are the ele- 
ments of the orthogonal transformation between the HO 
states (fi?72|(ni^i,n^)A) and (Ci^KA/'-C, 7V2L2)A). 



Combining the expansions of Eqs. (46) and (A3 1 it is 
possible to write the following expression for the A = 3 
"direct" potential kernel in the SNP basis: 
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(A4) 



r 



Finally, for the A = 3 system, the "exchange" part of of the interaction between the second-to-last and next- 

the potential kernel resembles closely the exchange part to-last nucleons (the target nucleons in this case), the 

of the norm kernel, and can be derived in a very similar matrix elements of which can be easily calculated using 

way as the latter. Indeed, besides different multiplicative the basis (44). Therefore, A — 3 "exchange" potential in 



factors, Eqs. (37) and (41) differ only for the presence the SNP basis is given by 



1 1 



1 i Si 

2 2 *i 



i[ Z s' 



[-2 T Til \ \ Z s\] y J a A 



z s' I ^ r 

\ I'l s'l J AAiCi I 



Ci Z s ^ \ Ci Z s 
J £ A j \ ^ /i si 

:{MiCiSihTi\V{V2^iaia2TiT2)\niiiSihTi) . 



{n'£',n[£[,A\AfiCi,ne,A)3 



(A5) 



Note that the above expression can be easily reduced 
to the exchange part of the norm kernel by replacing 
V{^/2£^iaia2TiT2) with 1. 



2. A> 4 

The expression derived in this Appendix are valid for 
systems with A>4: (a—l). 

We start by deriving the simplest of the integral ker- 
nels, i.e. the exchange part of the norm kernel (37). To 
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this aim, it is convenient to expand the (A — l)-nucleon 
eigenstates \A— 1 ail^^Ti) onto a HO basis containing 
anti-symmetric sub-clusters of A — 2 nucleons, e.g. 

\{NA^2'iA~2JA-2TA~2;nA~liA~ljA-l)hTi) . (A6) 

Here, the anti-symmetric states \NA-2iA-2JA-2TA-2) 
depend on the first A — 3 Jacobi coordinates of Eq. (13) 
(Cii '^2, • • • 7 Ca-s) and the first A— 2 spin and isospin co- 
ordinates, and are characterized by total number of HO 
excitations, spin, isospin and additional quantum num- 



bers Na-2, Ja-2, Ta-2, and iA-2, respectively. The ba- 



sis states (A6I are not anti-symmetrized with respect to 



the next-to-last nucleon, which is represented by the HO 

state {^A-2<:^A-irA-i\nA-i^A-ijA-i), where nA-i,£A-i 
are the HO quantum numbers corresponding to the har- 
monic oscillator associated with S,a-2, and jA-i is the 
angular momentum of the (yl — l)th nucleon relative to 
the cm. of the first A— 2. n terms of the basis states ( A6 1, 



the HO Jacobi channel state of Eq. (22 1 for the (^—1, 1) 
system can be written as 



l<i> 



= ^ {{NA~2iA^2JA-2TA-2; nA-liA-ljA-l)IlTi\A~l ai/f^Ti) 
X I [{{NA-^2iA-^2JA^2TA~2; nA^l^A^lj A-l)hTi \ ^^)«^; "-^] J^'^' 



(A7) 



where (^{NA-2iA-2JA-2TA-2',nA-i(-A-ijA-i)IiTi\A—\aiI'^^Tij are the coefficients of the expansion |20j 
of the (A— l)-cluster eigenstates on the basis (A6I, and the sum runs over the quantum numbers 

Na-2,t-a-2, Ja-2, Ta-2, nA-i,(-A-i, and j^-i- 

I 



According to Eq. ( 37 1 , in order to obtain the exchange 



part of the norm kernel we need to evaluate matrix ele- 
ments of the permutation corresponding to the exchange 
of the last two particles, Pa-i,a- The task can be accom- 



plished by, e.g., switching to a more convenient coupling 
of the nucleon quantum numbers (for a definition of the 
12-j symbol see Appendix Ib|) : 



[{{NA-2-iA-2JA-2TA-2;nA-iiA-ijA-i)IiTi; ^^)^T; n£] J'^T 
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. Ja-I (-A-l i S 



[NA-2iA-2JA-2\ {{7lA-l^A-lM)^S2)K]r) ( (T^_2 ^) ^ ) T) , 



(A8) 



and observing that, as a result of the action of Pa-i.a, 
the HO state {£,A-2'nA-i\{nA-i^A-i,nt}K) is changed 
into {S^A_2VA-i\{^A-i(-A-i,nt}K). The new set of Ja- 
cobi coordinates £,a-2 and ^a-i (obtained from i,A-2 and 
fjA-i, respectively, by exchanging the single- nucleon in- 



dexes A—1 and A) can be expressed as an orthogonal 
transformation of the unprimed ones. Consequently, the 
HO states depending on them are related by the orthog- 
onal transformation: 



{iA-2VA-i\i.nA-i^A-i,nC)K) 



^ {NL,MA-lCA-l,J^\nA-llA-l,nt,h.)A(A-2) 
NL,Ma-iCa-i 

x{-)L+^^-^-^{U-2ffA-imA-iCA-i,NL)A), 



(A9) 



26 



where the elements of the transformation are the general 
HO brackets for two particles with mass ratio d = A{A—2). 
After taking care of the action of Pa-i,a also on the spin 



and isospin coordinates, one can complete the derivation 
and write the following expression for the A > 4 exchange 
part of the norm kernel in the SNP basis: 
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r 



(AlO) 



where the second sum runs over the quantum numbers 

A'^A-2, *A-2, Ja~2, Ta-2, "A-I, ^A-1i JA-H ^A-l, 

and jA-i- The above expression was obtained expand- 
ing the 12-j symbol of Eq. ( A8 1 according to Eq. (Bl I or 



( B2 1 , and summing over the quantum numbers S2 and 



K2. Note that the norm kernel is symmetric under ex- 
change of primed and unprimed indexes and coordinates. 

We turn now to the derivation of the "direct" po- 
tential kernel of Eq. (40 1. As shown in Eq. pOl), in 



this case one needs to evaluate matrix elements of the 
interaction between the last two nucleons, (V^(fvi-i — 

rA,crA-iTATA-iTA){^ — Pa-i,a))- It is therefore useful 
to introduce two new Jacobi coordinates 



Ca-2 

Ca-1 



V2 



2 i^A-i + rA) - rA-2 



jA-i -rA) , 



(All) 
(A12) 



and switch to the HO basis states in which nucleons 
A—1 and A are coupled together to form two-particle 
states of the form ((a-iO'a-i<Jata-ita\N2L2S2J2T2), 
where N2 , L2 are the HO quantum numbers correspond- 
ing to the harmonic oscillator associated with Ca-1j and 
52, J2, and T2 are the two-nucleon spin, total angular 
momentum and isospin quantum numbers of the (A-1,A)- 
nucleons couple, respectively. This task can be achieved. 



e.g., continuing from the expansion of Eq. (A8I: 



Cl • • • £,A-2VA-lCri ■ ■ ■ CTA-lCTA 
x('?l • • • Ca-2Ca-10'1 • • • O-A-lf^ATl ■ ■ ■ Ta-iTa 



[NA-2^A^2JA-2;{{nA-liA-lM)^S2)K] r) (t^ ■ ■ ■ TA-iTa\({Ta-2\)Ti]^T 



'^f,kY^f2j2\^'f'^'] E 

T2,J2 I 2 )N2L2j 



^ r £ L2 A 1 

(-)^^+^ <^ \ {N2L2,AfC,K\n(., ua-uIa-uA) _^ 

Arc [ S2 K J2 ) 



[Na-2^A-2Ja-2Ta-2; {MC- N2L2S2J2T2)KT2] J^T 
I 



(A13) 



At this point, the expression for the "direct" potential kernel can be easily derived combining Eqs. (A6I and 



(A13l, and observing that Va-i,a{^ — Pa-i.a) is diago- 
nal in the quantum numbers NA-2'iA-2JA-2TA-2, J^C, 
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and S2J2T2 (for a definition of the 12-j symbol see Ap- 
pendix [bJ : 



(v4 - 1) ^ Rn'i'{r')Rne{r) ^ {{NA-2tA-2JA^2TA-2; n'A_i£'A_ifA-i)I[T[\A-la[l[''''T[) 

11' n 

x{{NA-2iA-2JA-2TA~2\nA-l^A~l3A~l)hTi\~l Qfi/f'Ti) 



xH 



l+2.]+I[+h+t'+l 



S2,J2.T2 I 2 J I 2 J 



2 ^ ^i: y V 2 
52 K Ja-2 1 r 5 ^2 if Ja-2 

J^i^^^A'^A^I 1 A' J i A J /i 



K A', A 



E E E<^2^2:-^^,A|nY,n:4_i,^:,_i,A')^(7V2L2,AA£,AKn^_i,^^_i,A)^ 



£ L'2 A' 1 f £ L2 A 



S2 K -h) I 5*2 K J2 



(A14) 



where the summation runs over the quantum numbers 
A''a-2, M-2, Ja-2, Ta-2, riA-i, (a-i, jA-i, as well as 
over the corresponding primed indexes. 

Finally we discuss the derivation of the "exchange" - 
potential kernel (41 1. The latter is a function of the 
matrix elements on the Jacobi channel states (A7l of 



the product of the P^-i.A exchange operator and the 
interaction between the {A — 2)th and {A — l)th nu- 
cleons: {^'l^,'^] pA-i,AVA-2,A-i\^in'^)- Therefore one 
may proceed, e.g., by first evaluating the action of Pa-i,a 
on the bra and then the matrix elements of 



Va-2,a-i between the modified bra and the ket j^;^^"^). 
For the first step one can utilize, as for the "exchange" - 
norm kernel, Eq. (A8). However, here, after the calcula- 



tion of the action of the exchange operator, it is conve- 
nient to perform the inverse of the transformation ( A8 1 
to return to the original coupling scheme of Eq. ( |A7| . In- 
deed, the interaction Va-2,a-\ acts on the (A-l)-cluster 
states and is diagonal in the quantum numbers n,£,s. 
The intermediate results resemble closely the expression 
of the exchange norm kernel and read: 
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1 rji/ rp '\ 



n'n 



{{NA^2iA-2JA^2TA-2\nA-lU-l]A~l)hTi\A-laJl'T^)f[f^{-f+^'^+'^' 




X JA-lJA-li'Jls's{-) 

t^A-\C.A-\JA-\ 



Ja-1 Z Ii 



/i Z s 



A 

La-i Z s' ) {nt,UA-iCA-iA\n'A-ii'A-i^n'e,K)A(A-2) 
i s J 

'^{{N'A_2i'A~2J'A-2T'A-2'MA-lJ^A-lJA-l)IlTi\VA~2,A~l\NA-2iA-2JA-2TA~2] nA-liA-liA-l)hTi), 

(A15) 




where the summation runs over both the primed and 



unprimcd sets of quantum numbers N'j: 



A-2' 'A-2' '^A-2' 



J'. 



"^'a-i^ ^a~i^ J'a-h and A'^a-2, M-2, Ja-2, Ta-2, 
n^-i, ^A-i, 2A-\- Note that, by replacing Va-2,a-\ 
with 1, one correctly recovers the exchange part of the 
norm kernel (AlO). For the second step, i.e. the eval- 
uation of the matrix elements of the interaction be- 
tween the second- and next-to-last nucleons, V(rA-2 — 
ta-i, crA-2<^A-iTA-2'''A-i)i we introduce two new Jacobi 
coordinates, namely 



PA-3- 



2{A-3) 



A-3 



A-1 



i=l 



AA-2 = (rvi-2 - 1 ) , 



(A17) 



and switch to the HO basis states in which nucleons 
A— 2 and A—1 are coupled together to form two-particle 
states of the form (pA-20'A-20■A-lTA-2TA-l|^^2^2S2i2^2), 
where n2 , ^2 s^re the HO quantum numbers correspond- 
ing to the harmonic oscillator associated with pA-2, and 
S2, J2, and ^2 are the two-nucleon spin, total angular mo- 
mentum and isospin quantum numbers, respectively: 



(^1 • ■ ■ Ca-3Ca-20'i • • • <J A-S.<^ A-2T1 ■ ■ ■ TA-3TA-2\{NA-2'iA-2JA-2TA-2', nA-ltA-lj A-l)IlTi) 
= '^{NA-ziA-zJA-^TA-Z; nA-2^A-23A-2\\NA~2iA-2JA-2TA-2) {-)'^^-^+'^^ fA-2 

x{-y^--+^^-' + J^-'+''JA-23A-2]A-lY.^ \ 51 ■^2j2t2(-p+*M 

Y I jA-1 II Y ) S2j2t2 I ^ Ti t2 

^A-2 h jA~2 



X 
A 



'Y^^ \ ^A^i \ jA-i\ X! I I (7^2^2,?^l^l, A|n^_i^A-i, "^71-2^^-2, A) 



S2 y j2 



\ S2 Y 

6 • • ■ PA-3PA-20-1 ■ ■ ■ (TA-20-A-1T1 ■ ■ ■ TA~2TA~l\{NA~3iA-3JA-3TA^3; ("-l^b n2(.2S232t2)Yt2) hTi) . (Al8) 

In deriving the above expression, we have expanded the (A-2)-nucleon anti-symmetric states \NA-2iA-2JA-2TA-2) 
onto a basis containing anti-symmetric sub-cluster of A — 3 nucleons, using the coefficient of fractional parentage 
(A^^_3i^_3 J^_3T/t_3; nyi_2^A-2jyi-2||A'A-2*yi-2'^A-27A_2)- The summation is intended over the quantum numbers 

Na-3, *A-3, Ja-3, Ta-3, nA-2, -^^-2, and jA-2- 
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In this basis, which is not anti-symmetric for ex- 
changes of the {A — 2)th nucleon, the anti-symmetric 
|A''^_3i^_3 J^_3r4_3) states depend on the first A — i 

Jacobi coordinates of Eq. (13 1 {^1,^2,- t^a-a) and 
the first A — 3 spin and isospin coordinates. Here 
Na-3, Ja-3, Ta-3, and iA-3 are total number of HO ex- 
citations, spin, isospin and additional quantum number 
characterizing the (A — 3)-nucleon anti-symmetric basis 
states, respectively. The second-to-last nucleon is repre- 
sented by the HO state {£,A-3<^A-2TA-2\nA-2iA-2jA-2), 
where nA-2T iA-2 are the HO quantum numbers corre- 
sponding to the harmonic oscillator associated with ^a-3i 
while jA-2 is the angular momentum of the (A— 2)th nu- 



cleon relative to the cm. of the first A— 3 nucleons. The 



summation in Eq. ( A18 1 runs over the quantum numbers 

Na-3, Ja-3, Ta~3, 'iA-3, nA-2, iA-2, and jA-2- Fur- 
ther, {n2i2,niii,X\nA-2iA-2,nA-l£A-l,>^)(A--l)/(A-3) 

are the general HO brackets for two particles with 
mass ratio d = {A — 1)/ {A — 3), which are the 
elements of the orthogonal transformation between 
the HO states (Ca-36i-2|("-a-2^a-2, J^a-i^a-^A) and 
{pA-3PA-2 I {niii , n2£2) A) . 

It is now trivial to complete the derivation of the 
"exchange" -potential kernel by complementing Eq. ( A15 1 
with the following expression: 



{{NA-2i'A^2J'A-2TA^2i-^A-lCA^lJA~l)hTi\VA-2,A-l\NA-2iA-2JA-2TA-2; nA-liA-ljA-l)hTi) 
= y^(^A-3U-3JA-3^A-3;"-A-2^A-2JA-2ll^A-2^A-2'^A-2'^A-2) 

x{NA^3iA-3JA-3TA-3;nA-2eA-2jA-2\\NA-2lA-2JA-2TA-2){-y'^-'+'^-'+^^-'+'^-' 



^i'A-2j A-2JA~lj A-iJ'a-2J A~2T'a^2Ta~2 ^ 

S2i2*2 



52^2^2 



rr\ 1 rpf \ C T' 1 

J^A-3 2 ^A-2 I I ^A~3 9 ^A- 



h t2 



h ^1 ^2 



A', A 



Ja-3 j'A-2 J'a-2 \ J JA-3jA-2 J A- 2 
Ja-i h Y iXjA-l h Y 




X ^ ^ {n'2i'2,nili,'X!\NA-lCA-l,n'A_2(-'A~2i^') 4^ ("2^2, A|rtA-l^A-l, »-A-2^A-2, A) 

rii^i n2^2 "2'?2 



\ h 4 A' ' 




f ^1 ^2 A ] 


152^^2 J 




\S2Y32\ 



{n'2^'2S232t2\V{^/2pA-2(^A-2f^A-lTA-2TA-l)\n2i2S2j2t2), 



(A19) 



were the summation runs over he quantum numbers 

Na-3-, *A-3, Ja-3-, Ta-3, nA-2, f.A-2, jA-2, "A-2' ^'a-2i 

and j'A_2- 

As for A ~ A {a = \) the {A — 2)-nucleon states 
\NA-2iA-2JA-2TA-2) are simply antisymmetric two- 
nucleon states of the kind \N2L2S2J2T2) characterized by 



a single Jacobi coordinate (^1), the transformation (AI8I 



is somewhat different for the four-nucleon system, leading 
to an independent expression for the matrix elements of 
the V2.3 interaction term between the target basis states 
(for a definition of the 12-j symbol see Appendix [B| : 
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((iV^45^J^r^;AA3£3J3)/iri|F2^3|A^2i2^2J2r2;n3%3)/iri) 

1 



l-(-l)^ 



1 1 ji/ ~1 (11 



To 



J^f^J^A'^A^ £3 J/ ^ 



r A', A 



2 2 «2 



^' -^2 'S'2 \ 



is J2 \ Si 
\ L2 82^ 



X (^2^2, "-1^1, A|n34, N2L2, A)3 



^1 4 A' 1 J 4 4 A 

S2 5^ j2 J I S2 J2 



^ E E E("24,m^i,A'|AA3/:3,A^24,A') 

ni£i n2^2 "2^2 



(»^2^2*2j2i2|^(\/2p2 Cr2 (73 T-2 T3) |n2^2S2j2i2) 



(A20) 



Note that in order to recover the full expression for the 
A — 4: {a—1) "exchange" -potential kernel, it is sufficient 
to replace A with 4 in Eq. (A15), and combine the latter 
equation with Eq. ( A20 1 . 



APPENDIX B: 12-j SYMBOL DEFINITION 



The 12-j symbol of the first kind [33] is defined by 



r- = « I = S^(^\a.+b+c+d+e+f+g+h+p+q+r+s-X'0-2\ ahX\\cdX\\efX\\ghX\ , , 
. \ a d \ ^ Icdp n«/9 lis'"- II ' 



(B2) 
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